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AFIT/GE/ENG/02M-09 


ABSTRACT 

This  research  presents  an  ambiguity  function  analysis  of  the  digital  audio 
broadcast  (DAB)  waveform  and  one  signal  detection  approach  based  on  signal  space  projection 
techniques  that  effectively  filters  the  direct  path  signal  from  the  receiver  target  channel. 
Currently,  most  Passive  Coherent  Eocation  (PCE)  research  efforts  are  focused  and  based  on 
frequency  modulated  (PM)  radio  broadcasts  and  analog  television  (TV)  waveforms.  One  active 
area  of  PCE  research  includes  the  search  for  new  waveforms  of  opportunity  that  can  be  exploited 
for  PCE  applications.  As  considered  for  this  research,  one  possible  waveform  of  opportunity  is 
the  European  digital  radio  standard  called  DAB. 

Por  this  research,  the  DAB  performance  is  analyzed  for  application  as  a  PCE  waveform 
of  opportunity.  Por  this  analysis,  DAB  ambiguity  function  calculations  and  ambiguity  surface 
plots  are  created  and  evaluated.  Signal  detection  capability,  to  include  characterization  of  time- 
delay  and  Doppler-shift  measurement  accuracy  and  resolution,  is  investigated  and  determined  to 
be  quite  acceptable  for  the  DAB  waveform. 

In  practice,  matched  filter  detection  may  be  used  in  a  PCE  receiver  for  exploiting  the 
DAB  waveform.  Por  this  signal-processing  scheme,  two  receiver  channels  are  required.  One 
channel  is  dedicated  to  establishing  a  good  reference  signal  as  derived  from  a  line-of-sight  signal 
reception.  The  second  channel  is  dedicated  to  target  detection  and  estimation,  i.e.,  the  antenna 
pattern  of  this  channel  is  directed  towards  a  specific  surveillance  region  and  receives  reference 
waveform  reflections  from  targets  of  interest.  The  correlation  response  between  the  reference 
waveform  and  the  target  return  (having  specific  time-delay  and  frequency  shift  characteristics) 
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can  be  used  with  an  appropriate  thresholding  seheme  to  deelare  target  presenee  and  estimate 
target  parameters. 

One  serious  problem  is  eneountered  with  the  proposed  teehnique,  namely,  the  existenee 
of  a  referenee  waveform  response  in  the  reeeiver  target  ehannel.  Without  eompensation,  this 
“direet  path”  eomponent  renders  target  deteetion  via  simple  matehed  filtering  nearly  impossible. 
However,  if  the  direet  path  eomponent  ean  be  eompletely  removed,  or  signifieantly  attenuated, 
matehed  filtering  ean  yet  be  a  viable  alternative  for  signal  proeessing.  This  researeh  proposes  a 
unique  time-domain  filtering  approaeh  to  eliminate  (or  signifieantly  reduee)  the  direet  path  signal 
eomponent  from  the  reeeiver  target  ehannel  of  a  PCL  system. 

The  speetral  estimate  of  a  reeeived  signal  eonsisting  of  a  direet  path  and  noise  eomponent 
is  first  determined.  In  this  ease,  the  direet  path  eomponent  eonsists  of  signals  from  two  different 
DAB  transmitters  in  a  Single  Frequeney  Network  (SFN).  The  reeeived  signal  is  then  filtered 
using  a  projeetion  operation  that  effeetively  projeets  the  reeeived  signal  into  the  referenee  signal 
null  spaee,  i.e.,  a  signal  spaee  ideally  eontaining  no  referenee  signal  energy.  Therefore,  the  only 
remaining  eomponents  of  the  reeeived  signal  response  lying  in  the  null  spaee  are  due  to  targets 
and  noise.  These  eomponents  are  then  eombined  and  a  post-filtered  speetral  estimate  formed. 
As  observed,  the  amount  of  referenee  signal  power  passing  the  filter  is  sharply  redueed.  This 
reduetion  in  referenee  signal  power  is  the  phenomenon  that  enables  the  target  response  of  the 
reeeived  signal  to  be  deteeted  using  a  post-filtered  detection  surface. 
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AMBIGUITY  FUNCTION  ANALYSIS  AND  DIRECT-PATH 
SIGNAL  FILTERING  OF  THE  DIGITAL  AUDIO  BROADCAST 
(DAB)  WAVEFORM  FOR  PASSIVE  COHERENT  LOCATION 

(PCL) 

CHAPTER  1 
INTRODUCTION 

I .  I  Background 

In  the  radar  world,  the  first  task  is  to  deteet  an  objeet  in  a  region  of  interest.  In  this 
eontext  the  objeet  to  be  deteeted  is  ealled  a  target.  Following  deteetion,  we  would  like  to 
determine  speeifie  target  parameters,  sueh  as  position  and  veloeity  and  if  possible,  we  would  like 
to  traek  the  target  and  its  related  parameters.  As  a  more  ehallenging  task,  we  would  also  like  to 
know  what  the  target  is,  i.e.,  identify  it  as  friend  or  foe. 

Deployed  radar  systems  employ  many  different  teehniques  to  aehieve  these  goals.  Many 
send  a  signal  speeifieally  designed  to  enhanee  the  signal  proeessing  tasks  of  deteeting,  traeking 
and  possibly  identifying  the  target.  The  radar  systems  of  eoneern  for  this  work  rely  on  target 
radio  frequeney  (RF)  re -radiation,  or  seattering,  eharaeteristies  to  aoeomplish  these  tasks.  For 
passive  eoherent  loeation  (PCF),  the  original  souree  of  energy  is  derived  from  a  waveform  of 
opportunity  in  the  environment  and  the  re-radiated  energy  from  the  target  is  proeessed. 

Both  aetive  and  passive  target  deteetion  teehniques  have  inherent  problems  assoeiated 
with  them.  Passive  teehniques,  whieh  rely  on  target  RF  emissions  sueh  as  target  radar  or  IFF 
transmissions,  ean  beeome  ineffeetive  in  deteeting  targets  whieh  implement  “emission  eontrol” 
(EMCON)  proeedures.  On  the  other  hand,  transmission  of  the  speeial  waveforms  sueh  as  those 
used  in  aetive  radar  systems,  may  be  intereepted  by  the  “Eleetronie  Counter  Measures”  (ECM) 
equipment  aboard  the  platform.  This  means  the  target  ean  detect  the  radar  before  the  radar 
deteets  the  target. 
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Another  problem  with  aetive  radar  systems  is  the  need  for  alloeated  frequeney  speetrum 
for  transmission  of  the  radar  waveform.  The  frequeney  speetrum  is  already  erowded  with  many 
eommunieation  and  broadeast  applieations.  Radar  waveform  therefore,  must  be  transmitted 
within  a  limited  frequeney  speetrum  as  dietated  by  other  applieations. 

As  a  ease  in  point,  radars  operating  at  “Very  High  Frequeney”  (VHF)  are  effeetive  in 
deteeting  hidden  ground  targets  sueh  as  tanks  in  a  forest.  These  radars  also  have  a  very  high 
range  eapability,  whieh  ean  be  important  for  deteeting  ballistie  targets.  However,  FM  radios  also 
broadeast  at  VHF  and  ean  affeet  the  performanee  of  the  radar  signifieantly.  In  the  end,  the 
benefits  assoeiated  with  VHF  operation  beeome  questionable  in  a  region  filled  with  FM  radio 
transmissions  [1]. 

In  faet,  the  above  deseribed  aetive  and  passive  teehniques  are  not  the  only  available 
deteetion  methods.  In  1978,  IBM  seientists  eame  together  in  an  attempt  to  enter  the  radar 
market.  As  part  of  their  effort,  seientists  looked  for  a  means  to  enhanee  target  deteetion  and 
traeking.  There  were  several  choiees  eonsidered: 

1.  Monostatie  aetive  radar:  In  this  teehnique,  the  radar  transmits  and  reeeives  a 
waveform  speeialized  to  ease  target  deteetion  and  traeking  proeessing.  The  transmitter 
and  reeeiver  are  eolloeated  or  are  separated  with  a  relatively  small  distanee  in  eomparison 
to  expeeted  target  distanee.  Praetieally  speaking,  these  systems  are  vulnerable  to  “Anti 
Radiation  Missiles”  (ARMs).  Also,  in  the  monostatie  arehiteeture,  if  the  target  platform 
has  an  ECM  system  onboard,  it  is  highly  possible  that  the  target  platform  ean  deteet  the 
radar  signal  before  the  radar  deteets  the  target  platform.  In  this  ease,  the  target  platform 
may  avoid  being  deteeted. 
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2.  Bistatic  radar  [2];  This  technique  requires  one  transmission  from  a  distant 
transmitter  site.  So,  the  transmitter  and  reeeiver  are  not  eolloeated.  In  this  ease,  the 
transmitter  is  still  vulnerable  to  ARM  attack  and  the  target  platform  with  ECM  capability 
can  still  detect  the  presenee  of  a  radar  and  possibly  avoid  it. 

3.  Non-eooperative  bistatic  radar;  This  system  uses  transmissions  from  another  radar 
system,  whieh  may  be  either  friendly  or  foe.  In  this  ease,  the  performanee  will  be  limited 
since  there  is  no  control  of  the  transmission.  The  advantage  of  this  technique  is  that  if  the 
original  radar  system  is  being  jammed,  this  system  will  be  minimally  affected. 

4.  Eleetronic  Support  Measures  (ESM):  ESM  are  passive  in  nature  and  rely  solely  on 
target  platform  transmissions  for  deteetion.  It  is  relatively  cheap  to  implement  and 
operates  covertly  since  no  transmissions  are  generated.  However,  ESM  performance 
diminishes  when  the  target  platform  applies  EMCON. 

East  but  not  least,  there  are  currently  a  significant  number  of  RE  transmissions  present  in 
the  environment  designed  for  non-radar  related  purposes,  e.g.,  TV  and  EM  broadeasts  or  mobile 
communication  signals.  Such  signals  may  be  exploited  for  deteetion,  tracking  and  possible 
identification  of  targets  -  this  exploitation  falls  under  the  umbrella  of  PCE  teehniques.  It  was 
this  idea  that  was  adopted  by  IBM  radar  researehers.  This  eoncept,  called  Passive  Coherent 
Location  (PCL),  is  what  lead  to  the  development  of  the  Silent  Sentry®  system  by  Eoekheed 
Martin  and  inspired  the  designs  of  many  other  systems  around  the  world. 

As  indieated  above,  the  PCE  concept  was  introdueed  by  researehers  at  IBM.  There  are 
several  PCE  definitions  available  in  the  open  literature.  Eoekheed  Martin  researchers  define  it  as 
a  system,  whieh  uses  everyday  broadeast  signals,  sueh  as  those  for  television  and  radio,  to 
illuminate,  detect  and  track  targets  [3].  Another  PCE  definition  comes  from  NASA  scientists 
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who  work  with  Lockheed  Martin  eolleagues  in  an  effort  to  make  traeking  of  launehed  spaee 
vehieles  eheaper.  NASA  defines  a  PCL  system  as  one  that  eapitalizes  on  the  energy  emitted  by 
eommereial  broadeast  stations,  sueh  as  radio  or  television  networks.  “Using  a  ground-based 
system  of  antennas,  reeeivers  and  signal  proeessors,  PCL  traeks  this  energy  as  it  is  refieeted  off 
an  aseending  launeh  vehiele,  allowing  its  aseent  trajeetory  to  be  plotted.”[4].  Another  PCL 
definition  taken  from  the  Journal  of  Eleetronie  Defense,  explains  PCL  as  an  eleetronie  system 
that  traeks  eivilian  radio  and  TV  broadeast  signals,  deteeting  aireraft  by  analyzing  the  slight 
disturbanee  in  the  eommereial  wavelengths  eaused  by  their  flight  [5]. 

1 .2  Different  PCL  systems 

The  basie  idea  of  PCL  is  elear;  exploitation  of  existing  RF  emissions  from  broadeast  or 
eommunieation  signals  in  such  a  way  to  detect,  track  and  also  identify  objeets  called  targets. 
Different  approaehes  have  been  taken  by  researehers  in  determining  whieh  RF  emissions  should 
be  used  in  PCF  applieations.  Most  PCF  literature  is  based  on  FM  broadeasts  and  analog  TV 
signals  [6-10].  The  speeifie  ehoiee  of  a  “waveform  of  opportunity”  [11]  greatly  affeets  signal 
proeessing  of  the  system.  For  example,  if  one  uses  the  audio  or  video  earrier  of  an  analog  TV 
signal  [8],  then  a  narrowband  system  is  present  and  Doppler  analysis  ean  be  performed  along 
with  “Direetion  of  Arrival”  (DOA)  measurements.  Doppler  and  DOA  measurements  are  used  to 
estimate  the  range,  position  and  veloeity  of  a  target. 

One  ean  also  estimate  target  parameters  sueh  as  range,  position  and  veloeity  through 
“Time  Differenee  Of  Arrival”  (TDOA)  teehniques.  TDOA  teehniques  are  somewhat  easier  to 
implement  sinee  matehed  filtering  ean  be  used  to  extraet  TDOA  information  from  the  reeeived 
signal.  This  may  not  be  possible  with  an  analog  TV  signal  beeause  of  the  waveform  strueture. 
In  an  analog  TV  signal,  there  exists  a  synehronization  frame.  This  short  period  synehronization 
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frame  makes  TDOA  measurements  ambiguous  at  approximately  every  9.0  km  [9].  There  are 
other  available  waveforms,  sueh  as  an  FM  radio  broadeast  signal,  whieh  are  viable  for  TDOA 
proeessing.  Here  the  bandwidth  of  the  signal  greatly  affeets  measurement  aeeuraey  [12],  A 
range  resolution  of  2.0  km  has  been  aehieved  using  FM  broadeasts  [9].  Silent  Sentry®  system, 
the  first  off-the  shelf  eommereially  available  PCL  produet  in  the  world,  uses  FM  broadeasts  as 
the  waveform  of  opportunity  [3]. 

Available  waveforms  of  opportunity  are  not  limited  to  only  analog  TV  and  FM 
broadeasts.  Many  people  have  eell  phones,  whieh  means  a  large  number  of  mobile 
eommunieation  signals  are  present  in  the  environment.  Other  possibilities  inelude  digital  TV  or 
digital  radio  broadeast.  Many  eountries  have  deeided  to  use  digital  broadeasts  as  future  standard 
for  audio  and  video  transmission.  Although  eurrent  usage  is  not  very  widespread,  the  audienee  of 
digital  broadeasts  is  steadily  inereasing  and  Europe  is  leading  the  world  in  this  effort.  For 
example,  the  “Digital  Audio  Broadeast”  (DAB)  waveform,  is  a  well-establish  standard  for  digital 
audio  broadeasts  and  is  starting  to  get  popular.  In  Britain,  60%  of  the  population  has  aeeess  to 
DAB.  In  Germany,  the  rate  is  30%,  in  Franee  25%,  in  Belgium  86%,  in  Sweden  85%,  and  in  the 
Netherlands  45%.  In  other  eountries  sueh  as  Turkey,  digital  audio  broadeast  implementation  is 
in  the  planning  phase  [13]. 

Looking  at  the  growing  developments  in  the  digital  domain,  espeeially  in  digital  audio 
broadeasting,  one  is  lead  to  the  possibility  of  using  DAB  waveform  in  PCL  applieations.  An 
eneouraging  faet  is  that  the  waveform  standard,  whieh  ean  be  found  in  open  literature,  is  a 
widely  aeeepted  standard  in  many  eountries.  Lor  Europe,  the  “European  Teleeommunieations 
Standards  Institute”  (ETSI)  has  issued  a  standard  for  digital  audio  broadeasting.  The  standards 
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are  published  in  a  doeument  ealled  ETS  300  401  [14],  The  potential  widespread  use  of  DAB 
makes  it  a  waveform  to  be  eonsidered  for  possible  PCL  applieations. 

1 .3  Researeh  problem  statement 

The  standard  DAB  signal  with  a  bandwidth  of  approximately  1.5  MHz  is  viable  for  TDOA 
ealeulations  in  a  PCL  applieation.  In  TDOA  signal  proeessing,  matehed  Liters  ean  be  used  to 
extract  TDOA  information  from  the  received  signal.  Before  building  an  actual  system  that 
operates  with  DAB,  it  would  be  best  to  first  test  DAB  waveform  properties  in  a  passive  radar 
scenario.  One  powerful  analysis  tool  for  characterizing  the  performance  of  any  waveform  in  a 
radar  system  with  a  matched  filter  receiver  is  the  “Ambiguity  function”.  An  “Ambiguity 
Surface”,  i.e.,  a  two-dimensional  plot  of  the  ambiguity  function,  indicates  what  performance  can 
be  expected  from  the  waveform  in  question.  As  addressed  under  this  research,  waveform 
performance  characterization  includes  determining  the  achievable  accuracy  and  resolution  in 
time-delay  and  Doppler  shift  measurements  of  the  received  waveform. 

Making  TDOA  measurements  with  a  PCL  system  requires  a  minimum  of  two  receiver 
channels.  The  “reference  channel”  is  dedicated  to  receiving  and  generating  a  reference  signal  to 
be  used  in  matched  filtering.  The  reference  signal  is  the  direct  break-through  of  the  broadcast 
transmitter  arriving  at  the  receiver  along  a  line-of-sight  path.  The  second  channel,  the  “target 
channel”  is  dedicated  to  receiving  returns  of  the  transmitted  waveform  of  opportunity  that  have 
reflected  from  objects  (targets)  in  the  environment.  Matched  filtering  in  PCL  applications 
includes  cross-correlating  received  signals  in  both  channels. 

One  problem  with  this  matched  filtering  process  arises  from  the  fact  that  the  channels  are 
not  perfectly  isolated,  i.e.,  one  component  in  the  reference  channel  is  the  target  return. 
Additionally  there  is  a  leakage  from  the  reference  signal  component  into  the  target  channel.  The 
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effect  of  the  reference  signal  response  in  the  target  channel  is  more  severe  mainly  because  the 
reference  signal  arrives  at  the  receiver  with  much  higher  power  than  the  target  return.  The 
effects  of  this  leakage  have  been  addressed  in  literature  [11].  The  reference  signal  component  in 
the  target  channel  causes  false  detection  in  the  receiver,  unless  this  component  can  be  removed 
or  significantly  reduced. 

This  research  proposes  a  unique  time-domain  filter  to  remove  or  significantly  reduce  the 
reference  signal  component  in  the  target  channel.  The  analytical  tool  used  for  this  research  to 
evaluate  performance  is  the  calculation  of  the  ambiguity  function  and  the  creation  of  “detection 
surface”  derived  from  ambiguity  surface. 

1 .4  Assumptions 

There  are  several  important  assumptions  made  to  be  able  to  concentrate  on  the  problem 
this  research  is  addressing. 

First,  the  geometry  for  the  PCL  system  in  consideration  is  monostatic.  Thus,  the  baseline 
in  the  transmitter-target-receiver  geometry  is  assumed  to  be  significantly  shorter  than  the  target 
distance  both  from  the  transmitters  and  receiver.  There  are  two  transmitters,  one  near  and  one  far 
away  from  the  receiver  site.  For  theoretical  calculations  and  simulations  requiring  a  target  return 
signal,  only  one  target  return  signal  is  embedded  into  the  received  signal.  Thus,  it  is  a  one-target 
scenario. 

The  next  important  assumption  is  the  condition  that  the  reference  signal  components  in 
both  reference  and  target  channels  arrive  at  the  receiver  simultaneously.  This  is  required  to  be 
able  to  correctly  calculate  the  time  delay  of  a  target  return  in  the  target  channel  of  the  matched 
filter  receiver  for  passive  radar  applications. 
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The  noise  eomponent  is  sampled  from  a  Gaussian  distribution  adjusted  for  eorreet  SNR  at 
the  reeeiver  input.  This  is  not  a  striet  assumption,  meaning  that  the  distribution  of  noise  does  not 
affeet  the  eonduet  of  this  researeh. 

Further  assumptions  are  stated  where  they  are  applied  throughout  this  researeh. 

1.5  Seope 

This  researeh  deals  with  a  PCL  system  operating  with  Digital  Audio  Broadeast  (DAB) 
waveform.  Speeifieally,  the  reeeived  signal  is  proeessed  and  the  output  of  the  matehed  fdter 
analyzed.  The  performanee  of  DAB  as  a  waveform  of  opportunity  for  a  PCL  system  is  evaluated. 
A  time-domain  filter  is  designed  to  solve  a  praetieal  problem  of  a  real-time  DAB  PCL  reeeiver. 

1.6  Materials  and  equipment 

Simulations  were  ereated  in  MATLAB®  Version  6.0,  from  The  Mathworks  Ine.,  Natiek, 
MA.  The  simulations  are  implemented  on  a  laptop  with  a  650  MHz  Intel  Pentium  III  proeessor. 

1.7  Thesis  organization 

Chapter  2  looks  at  the  world  of  PCL.  Different  passive  radar  applieations  are  introdueed 
ineluding  a  system  exploiting  DAB  signal  as  a  waveform  of  opportunity.  Next,  basie  knowledge 
about  the  DAB  signal  and  its  modulation  seheme  are  presented.  Chapter  3  deals  with  ambiguity 
funetion  analysis  of  DAB  waveform  and  introduees  a  unique  time-domain  filter  for  direet-path 
signal  filtering  in  PCL  systems  exploiting  DAB.  In  Chapter  4,  the  effieieney  of  the  direet-path 
filter  is  evaluated  through  MATLAB®  simulations.  Then,  Chapter  5  ineludes  eonelusions  and 
reeommendations  for  future  work.  Finally,  a  full  list  of  aeronyms  and  a  copy  of  the  MATLAB® 
code  developed  to  implement  the  filter  design  is  placed  into  the  Appendix. 


13 


A  list  of  references  and  vita  are  also  provided  at  the  end. 
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CHAPTER  2 


BACKGROUND 


2.1  Introduction 

This  chapter  presents  background  material  on  Passive  Coherent  Loeation  (PCL)  systems. 
PCL  systems  use  advantages  gained  by  waveform  eharaeteristies  of  the  signal  they  exploit  to 
deteet  and  traek  targets.  For  example,  FM  radio  broadeasts  are  in  VHF  region  of  the  frequeney 
speetrum.  In  this  speetral  region,  radar  systems  are  eapable  of  operating  at  long  ranges,  whieh 
eould  mean  a  long-range  eapability  is  possible  for  PCF  systems  operating  at  VHF.  The  benefits 
and  possible  drawbaeks  of  operating  in  this  frequeney  speetrum  are  diseussed.  Next,  present 
PCF  systems  are  diseussed  followed  by  Digital  Audio  Broadeast  (DAB)  related  issues  from  the 
open  literature.  The  ehapter  ends  with  a  diseussion  about  DAB  waveform  eharaeteristies  whieh 
are  neeessary  for  understanding  how  a  PCF  system  may  operate  with  a  DAB  waveform. 

At  this  time,  no  open  literature  eould  be  found  on  filtering  the  referenee  signal 
eomponent  from  the  target  ehannel  signal  (a  more  detailed  problem  statement  ean  be  found  in 
ehapter  one),  a  major  issue  addressed  in  this  researeh. 

2.2  Analog  TV  and  FM  radio  broadeast  waveforms 

Analog  TV  signals  have  a  bandwidth  of  about  6.0  MHz  and  are  typieally  loeated  in  the 
several  hundred  MHz  region.  In  the  frequeney  speetrum  of  analog  TV  signals,  there  are  two 
earriers  present,  one  for  the  vision  the  other  for  sound.  There  is  another  distinet  frequeney  band 
for  eolor  information  [15].  Typieally  a  wideband  speetral  response  provides  good  range 
resolution.  However,  the  speetral  shape  of  TV  signals  negates  the  advantages  gained  from  a 
wideband  system.  For  example,  the  reeeiver  matehed  filter  output  will  be  high,  not  only  as  a 
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function  of  true  target  range,  but  also  at  loeations  eorresponding  to  large  eross-eorrelation 
between  the  earners.  In  addition,  the  unambiguous  range  is  in  the  order  of  9  km  beeause  of  a 
periodie  synehronization  strueture  in  the  analog  TV  signal  speetrum  [9]. 

Beeause  of  the  above  faets,  it  is  not  eonvenient  to  do  range  measurements  with  analog 
TV  signals.  One  attempt  by  Griffiths  and  Long  dating  baek  to  the  1980’s  yielded  relatively  poor 
PCL  results  sinee  their  system  relied  mainly  on  range  measurements  [10].  Subsequent  work  on 
analog  TV  PCL  relied  on  Doppler  and  bearing  measurements,  as  used  by  Howland  [7,8]. 

The  lower  frequeneies  used  in  FM  radio  broadeasts  provides  many  distinet  advantages. 
The  HF  and  VHF  speetral  regions  eould  be  defined  as  low  frequeneies.  The  use  of  low 
frequeneies  allows  designing  radar  sensors  whose  performanee,  in  terms  of  deteetion,  is  better 
than  the  performanee  of  elassieal  mierowave  sensors.  Diffraetion  effeets  at  these  frequeneies 
allow  the  deteetion  of  hidden  targets  through  forests  or  behind  smooth  objeets.  At  these 
frequeneies,  penetration  of  eleetromagnetie  waves  through  foliage  is  also  exploited  in  SAR 
imaging  to  deteet  masked  targets  or  ground  moving  targets  whieh  eould  not  be  deteeted  using 
eonventional  SAR  imagery  [1]. 

Furthermore,  at  these  frequeneies  resonanee  effeets  may  make  low  RCS  targets  more 
visible.  In  the  resonanee  band  targets  aet  like  dipoles  yielding  an  eleetrieal  eurrent  when 
illuminated  that  is  independent  of  RCS.  The  eleetrieal  eurrent  forms  an  eleetrieal  field  eausing 
an  enormous  inerease  in  RCS,  independent  of  size  and  aspeet  of  the  target.  In  addition,  low 
frequeney  radars  are  more  immune  to  Anti  Radiation  Missiles  (ARMs).  This  faet  is  due  to 
praetieal  aeouraey  problems  in  ealeulation  of  the  exaet  loeation  of  a  low  frequeney  radar  site. 
Another  advantage  of  operating  at  low  frequeneies  is  the  ability  to  deteet  over-the-horizon 
targets.  Sky  waves  refieeted  by  the  ionosphere  enable  this  [1]. 
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However,  low  frequeney  regions  like  HF  and  VHF  bands  are  usually  dedieated  to 
oommereial  use  and  eommunieation  or  navigation  applieations.  PCL  systems  designed  to 
operate  in  this  frequeney  region  may  take  advantage  of  these  available  waveforms  without 
needing  speeifie  alloeated  frequeneies  for  this  purpose.  Of  all  PCL  systems  eited  in  open 
literature,  the  Manatash  Ridge  Radar  from  the  University  of  Washington  and  the  Silent  Sentry® 
system  from  Loekheed  Martin  have  low  frequeney  eomponents  operating  with  FM  radio 
broadeast  waveforms  [9,3]. 

2.3  PCL  Systems 

This  seetion  is  dedieated  to  providing  information  for  different  PCL  systems  from  the 
open  literature.  These  systems  work  mainly  with  FM  and  analog  TV  signals  and  have  different 
eharaeteristies  and  signal  proeessing  techniques. 

2.3.1  Manatash  Ridge  Radar  (MRR) 

The  following  discussion  about  MRR  is  mainly  from  a  paper  written  by  John  D.  Sahr 
and  Frank  D.  Lind  from  the  University  of  Washington  [9].  In  their  paper,  a  novel  method  for 
radar  remote  sensing  of  the  upper  atmosphere  is  described.  In  this  method,  FM  radio  broadcasts 
around  100  MHz  are  used  in  a  PCL  scenario.  FM  broadcasts  have  high  average  power  and 
possess  good  ambiguity  surface  properties.  Range  resolution  obtained  from  the  system  is 
approximately  1-2  km  which  is  perfectly  suitable  for  a  system  set  up  to  detect  atmospheric 
irregularities  in  the  upper  e-region,  a  huge  natural  phenomenon  scattered  over  a  wide  range. 
Doppler  resolution  is  claimed  to  be  1.0  m/s. 

In  the  MRR  system,  placing  the  receiver  at  a  site  completely  shielded  by  a  mountain 
range  solves  the  problem  of  dealing  with  a  direct  path  signal  from  a  line-of-sight  transmitter. 
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The  basic  configuration  includes  the  FM  transmitter,  a  receiver  to  pick  up  the  target  response  and 
another  receiver  to  process  the  direct  path  signal  as  needed  to  perform  the  signal  processing. 
Range  and  Doppler  are  the  measured  parameters  and  processing  is  done  non-real-time. 

As  mentioned  the  system  is  designed  to  detect  atmospheric  occurrences,  but  the  system 
has  also  been  successful  in  detecting  aircraft.  Reportedly,  the  system  can  detect  aircraft  and 
meteor  trails  at  distances  up  to  250  km. 

2.3.2  Lockheed  Martin’s  Silent  Sentry®  PCL  System  [3] 

Information  about  Silent  Sentry®,  the  only  off-the-shelf  PCL  product,  was  mainly 
obtained  from  a  document  published  by  Lockheed  Martin  scientists.  Research  on  Silent  Sentry® 
dates  back  to  the  1980’s  when  IBM  entered  the  radar  business  and  introduced  a  passive 
technique  which  later  became  known  as  PCL.  IBM’s  radar  division  was  later  sold  to  Lockheed 
Martin  where  the  second  generation  of  Silent  Sentry®  is  currently  under  development. 

The  basic  system  configuration  consists  of  receivers  using  up  to  six  FM  broadcast  station 
signals.  Fixed  as  well  as  portable  receiver  installation  is  possible.  Fundamental  system 
components  include: 

1 .  Target  Antennas  for  receiving  returns  from  the  surveillance  region. 

2.  Reference  Antennas  for  receiving  the  direct  path  reference  signal  from  the  FM 
transmitter. 

3.  High  Dynamic  Range  Receivers  as  needed  to  receive  the  strong  direct  path,  as  well 
as,  weaker  target  returns  arriving  at  the  same  time. 

4.  A/D  Converters. 

5.  Signal  processor  and  displays. 

6.  High  Speed  Tape  System  as  needed  for  fast  information  storage. 
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The  system  is  able  to  effeetively  traek  targets,  even  when  only  one  illuminator 
illuminates  the  target,  in  whieh  ease  traeking  is  two-dimensional  and  not  very  aeeurate.  Traeking 
quality  inereasing  with  more  FM  transmitters  illuminating  the  target,  whieh  leads  to  three- 
dimensional  traeking  using  three  and  more  illuminators. 

Mission  planning  tools  are  ineluded  in  the  system  to  enable  deployment  or  performanee 
simulations  around  the  world.  Advertised  system  performanee  ineludes  a  deteetion  range  up  to 
220  km  and  a  target  eapaeity  of  200  flying  objeets,  ineluding  everything  from  aireraft  to  ballistic 
missiles.  The  system  is  under  continual  development.  One  future  capability  may  include  the 
incorporation  of  TV  signals  not  limited  by  the  analog  format. 

2.3.3  TV  Based  Bistatic  Radar  [7,8] 

The  initial  TV  based  PCL  approach  of  Griffiths  and  Long  failed  to  provide  effective 
performance.  Since  that  time,  there  have  been  other  attempts  to  use  Direct  Broadcast  Satellite 
(DBS)  television  transmissions  for  range  measurement  using  TDOA  techniques  [16].  One 
successful  technique  used  analog  TV  waveforms  for  target  detection  and  tracking  and  was 
introduced  by  Howland  in  the  late  1990’s. 

Instead  of  using  range  measurements,  Howland  relied  on  Doppler  and  bearing 
measurements.  He  used  a  technique  borrowed  from  the  world  of  sonar,  i.e.,  target  location  is 
determined  by  measuring  both  Doppler  and  bearing  of  a  radiating  source  [17].  His  algorithm 
assumes  linear,  constant  velocity  target  movement  during  the  time  interval  when  location  is 
estimated.  Under  this  assumption,  a  given  Doppler  and  bearing  history  of  a  target  uniquely 
defines  its  location. 

The  signal  processing  in  Howland’s  system  begins  with  spectral  analysis  of  received 
signal,  where  targets  are  resolved  in  Doppler.  This  is  followed  by  a  time  Constant  False  Alarm 
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Rate  (CFAR)  detection  scheme  for  detecting  targets.  This  provides  an  indication  of  the  target’s 
presence  and  the  respective  Doppler  measurement.  Target  Direction  of  Arrival  (DOA) 
measurements  are  made  through  interferometry  using  two  Yagi  antennas  at  the  receiver.  A  4-by- 
4,  16-element  planar  array  to  provide  improved  performance  is  under  development.  A  Kalman 
filter  is  used  to  associate  target  Doppler  and  DOA  measurements.  The  association  of  a  target’s 
Doppler  with  its  DOA  information  is  important  for  subsequent  target  parametric  estimation. 
Doppler  and  bearing  profiles  belonging  to  each  target,  as  provided  from  the  association  of 
measurements,  are  used  in  an  extended  Kalman  filter  to  estimate  Cartesian  coordinates  and 
velocity  of  each  target.  Filter  initialization  is  done  using  a  Genetic  Algorithm  and  a  Levenberg- 
Marquardt  optimizer. 

Given  maximum  detection  and  tracking  range  is  approximately  260  km,  the  problem 
Howland  encountered  in  his  experiments  is  based  on  the  fact  that,  even  target  detection  was 
almost  complete,  the  system  could  only  track  one  third  of  the  detections.  The  problem  is  mainly 
attributed  to  the  low  information  content  of  the  received  target  response. 

2.3.4  Other  PCL  Concepts 

Clearly,  most  PCL  related  work  has  been  done  in  the  analog  signals  field.  However,  there 
are  some  research  efforts  using  complex  structured  digital  signals  as  PCL  waveforms  of 
opportunity,  including  cellular  phone  signals  and  digital  audio  and  TV  broadcasts  [18,19]. 

In  open  literature,  one  can  find  some  resources  on  the  exploitation  of  digital  broadcast 
waveforms  such  as  the  Digital  Audio  Broadcast  (DAB)  waveform.  This  waveform,  with  its 
relatively  wide  bandwidth  of  about  1.5  MHz,  offers  high  resolution  in  range  measurements 
through  correlation  of  the  received  signal  with  a  reference,  as  done  in  a  matched  filtering 
process.  In  the  DAB  standard  a  “Coded  Orthogonal  Frequency  Division  Multiplexing” 
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(COFDM)  technique  is  used  as  a  multiplexing  scheme,  which  provides  excellent  multipath 
resistance.  Experiments  using  the  DAB  waveform  as  a  PCL  waveform  show  promising  results 
[18,19]. 

In  theoretical  and  practical  experiments  using  the  DAB  waveform,  it  has  been  found  that 
severe  interference  from  the  direct  signal  (a  signal  arriving  at  the  receiver  along  a  line-of-sight 
path)  and  clutter  affect  DAB  performance  as  a  radar  waveform.  For  a  successful  detection 
scheme  via  matched  filtering,  the  direct  signal  component  received  in  the  receiver  target  channel 
should  be  eliminated  or  significantly  attenuated.  The  direct  signal  should  be  reduced  at  least  50 
dB  for  a  target  return  100  dB  below  the  direct  signal  level  to  be  detected  via  matched  filtering  at 
the  receiver.  The  reduction  is  also  needed  to  ensure  a  reasonable  dynamic  range  of  the  receiver 
[19]. 

At  the  same  time,  there  is  a  strong  multipath  component  in  the  received  waveform  that  is 
mainly  attributable  to  multiple  transmitters  using  the  same  waveform  as  done  in  the  “Single 
Frequency  Network”  (SFN)  structure  supported  by  DAB  signals. 

There  has  been  some  effort  put  into  removing  the  direct  signal  component  from  the 
receiver  target  channel.  There  may  have  been  some  unofficial  successful  results;  but,  no 
approach  has  been  publicly  discussed  in  open  literature.  For  this  research,  a  time  domain  filtering 
technique  based  on  signal  projection  spaces  is  proposed  for  removing  the  direct  signal 
component  from  the  target  channel  of  the  PCF  receiver.  For  this  implementation,  it  is  important 
to  understand  the  basics  of  the  DAB  waveform  and  the  COFDM  technique.  The  next  section 
provides  relevant  information  on  the  important  aspects  of  those  two  concepts. 
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2.4  Digital  Audio  Broadcast  (DAB)  waveform 


The  digital  audio  technology  enabling  CD  quality  sound  was  first  developed  for  satellite 
delivery  in  the  early  1980s.  The  system  operated  in  10-12  GHz  frequeney  region  and  was  not 
suitable  for  mobile  reception.  In  addition,  “local  services”,  an  important  concept  of  FM 
broadeasting,  was  not  possible  via  satellite  delivery  and  terrestrial  digital  radio  broadeasting 
emerged  as  an  alternative.  After  deeiding  that  earth-based  transmission  was  the  best  approach, 
choosing  the  right  modulation  scheme  and  the  right  frequency  interval  beeame  an  issue.  Initial 
work  on  using  Pulse  Code  Modulation  (PCM)  in  the  existing  FM  broadcast  frequency  band 
failed  due  to  interference  and  spectral  inefficiency  [13]. 

After  years  of  researeh  and  development,  a  consortium  called  Eureka  147  established  a 
standard  for  terrestrial  digital  radio  broadeasting  called  Digital  Audio  Broadeasting  (DAB).  The 
chosen  DAB  transmission  method  was  Coded  Orthogonal  Frequency  Division 
Multiplex  (COFDM).  More  than  10  years  after  the  first  offieial  presentation  of  DAB  at  the 
World  Administrative  Radio  Conferenee  1988  (WARC’88),  DAB  is  still  in  the  introduetion 
phase.  There  are  some  countries  across  Europe  with  existing  DAB  transmitters  serving  the  public 
but  in  most  places  the  system  is  still  under  test  or  in  the  planning  phases  [13]. 

The  core  of  the  DAB  transmission  system  is  COEDM,  a  multi-carrier  signal  transmission 
technique  that  uses  multiple  low  data  rate  subearriers  to  transmit  high  rate  data  streams.  The 
next  section  introduces  the  COEDM  concepts  important  to  understanding  this  research. 

2.4.1  The  COEDM  Modulation  Teehnique  [20] 

The  COEDM  teehnique  is  a  speeial  applieation  of  “Orthogonal  Erequency  Division 
Multiplex”  (OEDM)  where  “Eorward  Error  Correction”  (EEC)  is  used  to  further  reduce  the  “Bit 
Error  Rate”  (BER)  in  a  digital  eommunieation  system.  OEDM  can  be  viewed  as  either  a 
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modulation  or  multiplexing  technique.  The  main  reason  for  using  OFDM  is  its  robustness 
against  frequency  selective  fading  or  narrowband  interference  resulting  from  its  multi-carrier 
structure.  A  single  faded  frequency  or  interference  source  at  a  particular  frequency  only  affects  a 
few  of  the  multi-carriers,  also  called  “subcarriers”  in  OFDM. 

There  are  different  ways  to  implement  “frequency  division  multiplexing”  (FDM).  The 
first  technique  involves  dividing  the  total  available  channel  bandwidth  to  N  non-overlapping 
frequency  sub-channels.  This  technique  leads  to  spectral  inefficiency;  therefore,  overlapping  of 
individual  sub-channels  is  proposed  to  provide  more  efficient  bandwidth  use  at  the  expense  of 
increasing  BER.  To  increase  system  spectral  efficiency  further,  while  at  the  same  time  reducing 
BER,  cross-talk  between  the  subcarriers  is  prevented  (minimized)  by  choosing  orthogonal 
subcarrier  frequencies  [20]. 
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The  frequency  bin  of  the  subcarrier  with  f^.  as  center  frequency 


Figure  2.1  Orthogonal  Subcarriers  in  OFDM 


If  the  OFDM  subcarriers  are  orthogonal,  as  shown  in  Figure  2.1,  the  receiver  can  be 
implemented  to  act  as  a  bank  of  demodulators.  In  this  case,  the  receiver  translates  each 
individual  subcarrier  down  to  baseband  and  integrates  the  resultant  signal  over  a  symbol  period 
to  recover  the  raw  data.  As  each  carrier  is  down-converted  to  baseband,  the  other  orthogonal 
subcarrier  responses  in  that  bank  will  have  an  integer  number  of  cycles  in  the  symbol  period  Ts, 
which  effectively  integrates  to  zero.  Thus,  the  subcarriers  are  linearly  independent  (orthogonal) 
if  the  carrier  spacing  is  an  integer  multiple  of  1/1$. 

The  DAB  signal  consists  of  a  sum  of  subcarriers  that  are  independently  modulated  using 
“Differential  Quadrature  Phase  Shift  Keying”  (DQPSK).  If  z/  are  the  complex  DQPSK  symbols. 
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K  is  the  number  of  OFDM  subearriers,  Ts  the  symbol  duration  and  fc  the  earrier  frequeney,  then 
one  OFDM  symbol  starting  at  t  ean  be  written  as: 

[kh-x  r  f  /  +  0  5^  11 

s{t)  =  Re]  ^  exp  jlTt  /, - —  (t  - t  +  7^ 

[i=-KI2  L  V  )  JJ  (2.1) 

s{t)  =  0  t  <tg  or  t>t^+Tg 

In  literature,  the  eomplex  baseband  signal  representation  is  used  as  given  in  Equation  2.2. 

KI2-\  f 

■^(0=  X  ^,+jf/2exp  j27U  —  {t-ts) 

i=-KI2Y_  V 

s{t)  =  0  tg>t>tg+Tg 

The  complex  baseband  OFDM  signal  representation,  as  given  in  Equation  2.2,  is  nothing 
more  than  the  inverse  Eourier  transform  of  K,  DQPSK  input  symbols.  The  time  discrete 
equivalent  is  the  inverse  discrete  Eourier  transform  (IDET)  as  given  by  Equation  2.3. 

f  .  .  in^ 

s(n)  =  ZiBxp  j2;r—  (2.3) 

I  NJ 

In  Equation  2.3,  N  is  the  total  number  points  used  in  the  IDET  process.  In  practice,  this 
transform  can  be  implemented  most  efficiently  using  an  “inverse  fast  Eourier  transform”  (lEET). 
To  do  this,  the  vector  of  complex  DQPSK  symbols  is  zero-padded.  Zero-padding  helps  to 
prevent  intolerable  aliasing  when  the  IFET  output  samples  are  subsequently  passed  through  a 
digital-to-analog  converter.  By  zero-padding,  the  vector  of  complex  DQPSK  symbols  is 
lengthened  to  a  power  of  two.  Eor  the  DAB  waveform  there  are  1536  subcarriers  and  after  zero¬ 
padding  the  complex  symbols  vector  has  a  length  of  2^'  =  2048. 
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Still,  the  IFFT  output  is  not  ready  to  transmit.  To  eliminate  “intersymbol  interference” 
(ISI)  in  the  system,  a  guard  time  is  introduced  into  each  OFDM  symbol.  Here,  a  symbol  is 
defined  as  the  IFFT  output  or  the  sum  of  all  DQPSK  data  modulated  subcarriers.  The  guard 
interval  is  chosen  to  be  larger  than  the  expected  delay  spread,  such  that  multipath  components 
from  one  symbol  cannot  interfere  with  the  next  symbol. 

The  guard  interval  could  consist  of  no  signal  at  all,  i.e.,  the  transmitter  stops  transmitting. 
In  this  case,  “inter-carrier  interference”  (ICI)  between  the  subcarriers,  also  called  cross-talk, 
could  occur  due  to  lost  orthogonality  between  subcarriers.  This  effect  is  illustrated  in  Figure  2.2. 
When  the  OFDM  receiver  tries  to  demodulate  the  first  subcarrier,  it  will  encounter  some 
interference  from  the  second  subcarrier  because,  within  the  FFT  interval,  there  is  no  longer  an 
integer  number  of  cyclic  difference  between  subcarrier  1  and  2. 


Part  of  subcarrier  #2  causing 
ICI  on  subcarrier  #I 


Subcarrier  #I 


OFDM  symbol  time 


Figure  2.2.  Effect  of  multipath  using  no  signal  during  the  guard  iterval 


26 


To  eliminate  (minimize)  ISI  without  eausing  ICI,  the  OFDM  symbol  is  oyolieally 
extended  during  the  guard  time.  This  proeedure,  ealled  eyelie  prefix  insertion,  is  implemented 
by  appending  a  eopy  of  the  last  A  seeonds  (or  equivalent  number  of  samples)  of  the  IFFT  output 
to  the  beginning  of  the  OFDM  symbol,  where  A  is  the  guard  time  duration.  Cyelic  prefix 
insertion  ensures  that  delayed  replieas  of  the  original  OFDM  symbol  always  have  an  integer 
number  of  cyeles  within  an  FFT  interval,  as  long  as  the  delay  is  smaller  than  the  guard  time. 
Therefore,  multipath  signals  with  delays  smaller  than  the  guard  time  eause  neither  ISI  nor  ICI 
under  ideal  eonditions. 


Figure  2.3.  OFDM  Symbol  with  Cyelie  Prefix  Insertion 
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The  importance  of  cyclic  prefix  insertion  and  the  guard  interval  becomes  more  evident  in 
looking  at  the  DAB  transmission  characteristics.  The  DAB  standard  is  built  around  a  Single 
Frequency  Network  (SFN)  architecture.  In  a  SFN  transmitters  in  a  given  region  transmit  the 
same  signal  over  a  specified  channel,  hence,  multiple,  time-delayed  copies  of  the  same  signal  are 
simultaneously  received  which  may  be  correctly  called  multipath.  These  multipath  signals  do 
not  cause  ISI  if  they  arrive  at  the  receiver  with  a  delay  no  longer  than  the  guard  interval. 
Therefore,  the  guard  interval  duration  has  an  immediate  influence  on  the  type  of  SFN  that  can  be 
supported.  The  longer  the  guard  interval,  the  greater  the  maximum  allowable  distance  between 
SFN  transmitters  becomes.  However,  the  signal  energy  transmitted  during  the  guard  interval 
(cyclic  prefix)  is  not  used  by  the  receiver  for  data  demodulation.  Therefore,  using  a  guard 
interval  actually  decreases  the  channel  capacity.  Typically,  the  guard  interval  is  chosen  to  be  no 
greater  than  !  A  [21].  For  DAB  Mode  I  [14],  the  symbol  duration  is  1.0  ms  and  the  guard 
interval  is  chosen  to  be  0.246  ms. 


2.4.2  The  DAB  standard 

The  DAB  signal  structure  was  standardized  by  ETSI.  Details  of  the  standard  can  be  found 
in  ETS  300  401,  a  document  published  by  ETSI  about  DAB  related  issues.  In  this  section,  the 
DAB  signal  structure  is  given  as  stated  in  ETS  300  401: 

5(0  =  Re]  (2.4) 

m  =  -oo  /  =  0  k  =  -KI2 


0  for  1  =  0 

^  .Rect(t/Ts)  for  I  =  1,2,...,L 

Tg  =Tu+ A 

Various  DAB  signal  parameters  and  variables  are  defined  as  follows: 


(2.5) 

(2.6) 
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L  the  number  of  OFDM  symbols  per  transmission  frame  (exeluding  the  Null 
symbol); 

K  the  number  of  transmitted  subearriers; 

Tp  the  transmission  frame  duration; 

'^NULL  the  Null  symbol  duration; 

Tg  the  duration  of  OFDM  symbol  of  indices  /  =  1,  2,  3,...,Z  ; 

T^j  the  inverse  of  the  subcarrier  frequency  spacing; 

A  the  duration  of  the  guard  interval; 

I  k  complex  DQPSK  symbol  associated  with  carrier  k  of  the  OFDM  symbol  I 
during  transmission  frame  m  .  Its  values  are  defined  in  the  following  sub-cases. 
For  A:  =  0,  ^  =  0,  so  that  the  central  carrier  is  not  transmitted; 

the  central  frequency  of  the  signal.  [14] 

The  DAB  Mode  I  transmission  frame,  as  used  in  this  research  analyses,  includes  a  total  of 
77  symbols.  The  first  symbol  is  the  “Null  Symbol”.  During  transmission  of  the  null  symbol,  the 
composite  signal  s{t)  equals  zero.  The  Null  symbol  and  the  following  “Phase  Reference 
Symbol”  together  represent  what  is  called  the  Synchronization  Symbol  and  are  used  by  the 
receivers  to  synchronize  to  each  transmission  frame.  Each  DAB  communication  symbol  consists 
of  1536  subcarriers  and  is  1.246  ms  long,  including  a  cyclic  prefix  of  0.246  ms.  Only  the  null 
symbol  has  a  duration  of  1.297  ms. 

The  phase  reference  symbol  is  not  only  used  as  part  of  the  synchronization  process  but 
also  as  a  reference  for  differential  data  modulation.  The  phase  reference  symbol  is  defined  by  the 
values  of  z,  ^  for  /  being  the  phase  reference  symbol  index. 
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for  — K/2<k<0  and  0  <  k  <  K  /  2 
for  k  =  0 


(2.7) 


^i,k 


0 


<Pk  =  (2.8) 

The  phase  referenee  symbol  is  ereated  at  the  reeeiver  using  a  known  sequenee  ,  whieh 
is  a  funetion  of  transmission  mode  and  earner  index.  Eaeh  is  an  integer  number  in  the 

interval  [0,3]  and  is  used  to  modulate  eaeh  subearrier.  Eaeh  number  eorresponds  to  a  two-bit 
pattern.  When  the  signal  arrives  at  the  reeeiver,  the  phase  value  for  eaeh  subearrier  is  found  using 
East  Eourier  Transform  (EET).  As  explained  previously,  the  Inverse  East  Eourier  Transform 
(lEET)  has  been  used  in  the  transmitter  to  ereate  the  transmission  waveform. 

Nevertheless,  the  reeeived  signal  is  not  a  funetion  of  a  unique  transmission  signal  from 
one  and  only  one  transmitter.  In  the  SEN  arehiteeture,  in  a  region  of  interest,  multiple 
transmitters  transmit  identieal  signals.  Thanks  to  the  guard  interval  in  the  signal  strueture 
different  transmitter  signals  at  the  reeeiver  appear  as  multipath  for  the  reeeiver  during  symbol 
time  Ts-  There  is  also  a  ehannel  noise  eomponent  in  the  reeeived  signal  whieh  will  effeet  the 
phase  determination  by  the  reeeiver  during  DAB  demodulation. 

To  make  the  situation  elearer,  a  ease  is  eonsidered  where  only  two  DAB  transmitters  are 
in  the  region.  Again,  the  transmitted  signal  in  a  given  ehannel  is  identieal  for  both  transmitters. 
At  the  reeeiver,  the  reeeived  signal  eonsists  of  a  strong  DAB  signal  eomponent  from  the  nearest 
transmitter,  a  weaker  DAB  signal  eomponent  from  the  more  distant  transmitter  and  a  noise 
eomponent.  The  signal  from  the  distant  transmitter  will  be  a  delayed  eopy  of  the  stronger  signal 
from  the  nearer  transmitter.  At  this  time,  it  is  useful  to  consider  the  complex  phasor 
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representation  of  individual  subearriers.  Speoifieally,  in  the  transmitter  every  subearrier  phase  is 
ehanged  aeeording  to  the  information  bits  and  each  subcarrier  may  be  represented  by  a  complex 
phasor  having  the  proper  phase  as  a  result  of  DQPSK  data  modulation.  In  the  receiver,  for  a 
given  signal  component,  the  received  phase  of  the  complex  phasor  for  each  subcarrier  is  a 
function  of  the  transmission  phase,  subcarrier  wavelength  and  propagation  distance  (delay  time) 
from  the  transmitter  to  the  receiver. 

In  demodulation  at  the  receiver,  the  signal  is  taken  to  baseband  and  Fourier  transformed, 
yielding  the  complex  amplitudes  for  each  subcarrier.  In  Figure  2.4  two  cases  are  indicated.  In 
case  a)  the  phasor  representation  for  individual  components  of  three  subcarriers  are  given  and  in 
case  b)  the  sum  phasor  as  a  result  of  the  FFT  at  the  receiver  is  given.  The  resultant  phasor  for  a 
given  subcarrier  after  FFT  is  the  vectoral  sum  of  each  component  shown  in  case  a). 


31 


i  Noise 


Sum  phasor  at  receiver 


i  Noise  i  Far  transmitter 

Near  transmitter 

t 

y 

Subcarrier  no  (  I  ) 

(ii) 

(III) 


Figure  2.4.  Phasor  representation  of  demodulated  subcarriers  (I-III) 


If  the  signal  from  the  nearest  transmitter,  which  is  the  strongest  component  in  the  receiver 
channel,  is  the  only  component,  data  demodulation  is  relatively  unaffected.  This  strongest  DAB 
component  in  the  receiver  from  the  nearest  transmitter  is  called  the  reference  signal.  In  a  given 
transmission  frame  for  this  reference  signal  the  phase  reference  symbol  is  demodulated  first.  As 
mentioned  earlier,  the  specific  two-bit  phase  pattern,  coded  into  each  specific  subcarrier  is 
known  by  the  receiver.  So,  after  the  FFT  is  applied  by  the  receiver,  the  resultant  phase  values  for 
each  subcarrier  are  matched  with  corresponding  two-bit  patterns  of  the  phase  reference  symbol. 
The  information  contained  in  each  transmitted  DAB  signal  includes  interleaved  time  and 
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frequency  data.  Therefore,  the  next  step  involves  deinterleaving  the  bits  in  time  and  frequency 
followed  by  channel  decoding  using  the  Viterbi  algorithm  [13]. 

In  reality,  the  received  signal  may  contain  a  strong  multipath  component  as  a  result  of  a 
different  DAB  transmitter  in  the  SFN.  In  this  case,  the  FFT  at  the  receiver  will  yield  a  complex 
phasor  representing  the  vectoral  sum  of  the  reference  signal  component  and  the  multipath 
component  for  each  subcarrier.  If  the  noise  component  is  also  taken  into  consideration,  a 
condition  similar  to  that  plotted  in  Figure  2.4  is  present.  For  the  receiver  to  perform  properly  the 
absolute  phase  value  of  each  subcarrier  phasor  is  unimportant.  Using  the  phase  reference  symbol 
the  receiver  creates  initial  phase  thresholds  for  demodulation.  It  makes  no  difference  if  the  phase 
thresholds  are  based  on  the  reference  symbol  phasor  only  (the  ideal  case)  or  on  the  sum  of  the 
reference  and  multipath  phasor,  as  long  as  the  SNR  remains  the  same.  As  seen  in  Figure  2.4,  the 
noise  component  also  affects  the  final  relative  phasor  phase  values  aj  ,a2  ,  and . 

In  Figure  2.5  an  example  is  shown  for  three  subcarriers.  In  case  a)  the  reference  signal 
phasor  and  multipath  component  phasor  are  added  vectorially,  yielding  the  sum  phasor.  The 
noise  component  is  shown  separately.  If  there  is  only  a  reference  signal  present,  the  SNR  is  a 
function  of  the  ratio  between  the  amplitudes  of  reference  signal  component  and  the  noise 
component.  However,  due  to  multipath  the  demodulation  will  be  a  function  of  the  final  sum 
phasor  given  in  Figure  2.5  and  the  SNR  becomes  a  function  of  the  amplitude  ratio  of  the  sum  of 
reference  and  multipath  phasor  to  the  noise  component. 


33 


i  Noise 


Multipath  component 


Reference  signal 


^  Sum  phasor  at  receiver 
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Figure  2.5.  Phasor  representation  of  demodulated  subcarriers  (I-III)  amplitude  comparison 


Case  b)  in  Figure  2.5  provides  a  comparison  between  the  reference  phasor  and  the  sum 
phasor  for  three  different  subcarriers.  For  subcarrier  (I)  and  (III),  the  sum  phasor  has  a  larger 
amplitude  and  yields  a  higher  SNR  resulting  in  improved  BER  before  FEC. 

SNR  improvement  due  to  addition  of  reference  and  multipath  components  at  the  receiver 
can  be  theoretically  analyzed.  An  example  is  given  in  Eigure  2.6.  In  Eigure  2.6,  the  reference 
and  multipath  phasors  are  represented.  The  amplitude  of  the  multipath  signal  is  smaller  than  the 
reference  signal.  The  phase  of  the  multipath  phasor,  relative  to  the  reference  phasor,  will  be 
uniformly  random  in  the  interval  [0,2;r].  Eor  a  given  multipath  amplitude  interval  [0,x]  the  tip 
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of  the  multipath  phasor  in  the  figure  must  fall  into  the  patterned  area  in  the  dashed  circle  with 
radius  x  to  result  in  a  sum  phasor  smaller  than  the  reference  phasor,  indicating  lower  SNR  in  the 
receiver.  The  patterned  area  region  in  the  dashed  circle  is  smaller  than  one-half  of  the  dashed 
circle.  Therefore,  it  is  more  likely  that  the  resultant  sum  phasor  has  a  higher  SNR  value. 


Reference  phasor 


Multipath  phasor 


Sum  phasor 


Figure  2.6.  Amplitude  comparison  of  sum  phasor  to  reference  phasor 


In  reality,  some  subcarrier  multipath  will  result  in  a  smaller  sum,  degrading  system  BER 
for  that  particular  subcarrier.  However,  in  a  larger  number  of  subcarriers  the  multipath  will 
result  in  a  larger  sum  and  improved  BER.  It  should  be  noted  that  the  system  also  includes  EEC 
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on  time  and  frequency  interleaved  bits.  Therefore,  it  can  be  concluded  that  multipath  signals  at 
the  receiver,  as  caused  by  SFN  transmission,  improves  overall  DAB  receiver  performance 
provided  the  multipath  delay  is  less  than  the  guard  interval  of  the  system. 

2.5  Summary 

In  this  chapter,  background  material  on  PCL  systems  is  presented.  Different  PCL 
systems  operating  with  different  waveforms  of  opportunity  are  introduced.  It  has  been  noted  that 
current  systems  mainly  exploit  analog  waveforms  such  as  FM  radio  broadcasts  and  analog  TV 
signals.  The  advantages  of  using  low  frequency  modulation  schemes  like  FM  are  stated. 

In  addition  to  analog  waveforms,  newer,  more  complex  waveforms  having  digital 
modulation  are  also  considered  for  PCL  applications.  Among  these,  the  Digital  Audio  Broadcast 
(DAB)  waveform  shows  promising  results.  A  closer  look  at  the  modulation  and  multiplexing 
scheme  of  the  DAB  architecture  is  shown.  It  has  been  found  and  shown  that  the  modulation  as 
well  as  the  multiplexing  scheme  does  support  possible  PCL  exploitation  of  the  waveform. 

In  the  next  chapter,  a  filter  design  is  proposed  to  eliminate  the  reference  signal 
component  in  the  PCL  receiver  target  channel.  In  addition,  the  ambiguity  function  is  explored  as 
a  means  of  analyzing  the  performance  of  the  DAB  signal  as  a  PCL  waveform.  Originating  from 
the  ambiguity  function,  an  analytical  tool  called  a  detection  surface  is  introduced  to  evaluate  the 
filter  performance. 
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CHAPTER  3 


METHODOLOGY 


3.1  Introduction 

In  this  chapter,  the  tools  required  for  analyzing  the  Digital  Audio  Broadeast  (DAB) 
waveform  for  PCL  applieations  are  introdueed.  The  diseussion  starts  with  a  thorough 
explanation  of  the  ambiguity  funetion,  a  powerful  analysis  tool  used  to  study  radar  waveform 
properties.  Next,  issues  related  to  the  DAB  waveform  ambiguity  funetion  are  stated. 
Afterwards,  a  time  domain  filter  design  is  introdueed  as  a  means  to  enable  time  differenee  of 
arrival  (TDOA)  analysis  of  the  DAB  signal  in  a  PCL  system.  The  ehapter  eoneludes  with  an 
explanation  of  the  signal  proeessing  steps  taken  to  implement  filtering  using  MATLAB®. 

3 .2  Ambiguity  Eunetion  Analysis  [  1 2] 

“The  ability  of  a  radar  to  deteet  targets  amongst  returns  from  other  objeets  (elutter)  in  a 
region  of  interest  and  the  ability  to  determine  parameters  of  these  targets,  sueh  as  range,  bearing, 
size  and  veloeity,  depends  largely  on  the  radar  waveform”  [11]. 

Eor  PCL  applieations  using  waveforms  of  opportunity,  one  does  not  have  the  ability  to 
synthesize  optimal  waveforms.  The  approaeh  is  to  use  eommereial  eommunieation  and  other 
broadeasted  waveforms,  already  present  in  the  environment,  as  radar  waveforms  for  deteeting, 
traeking  and  possibly  identifying  targets.  These  eommereial  signals  oeeupy  a  eertain  speetral 
region  and  have  distinet  modulation  oharaeteristies  assoeiated  with  them.  Unless  a  radar 
waveform  is  speoifieally  designed  to  work  for  oommunieation  and/or  broadeasting  purposes  one 
is  not  able  to  use  a  radar-optimized  signal  in  PCL  applieations.  In  determining  the  extent  to 
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which  a  waveform  of  opportunity  can  be  exploited  for  radar  implementation,  the  ambiguity 
funetion  and  resultant  ambiguity  surfaee  plot  provide  a  useful  tool  for  analysis. 

The  ambiguity  function  surface  plot  and  corresponding  calculations  can  be  used  to 
qualitatively  assess  ‘how  well’  a  partieular  waveform  aehieves  and  affects: 

1.  Target  detection, 

2.  Parametric  measurement  aeeuraey, 

3.  Range  and  Doppler  resolution, 

4.  Ambiguities,  and 

5.  Clutter  rejection. 

The  ambiguity  funetion  is  widely  used  as  a  tool  for  analyzing  radar  waveform 
performance  [12].  Basically,  it  is  a  series  of  correlation  integrals,  or  the  absolute  value  squared, 
based  on  matehed  filter  detection,  where  a  reeeived  signal  having  different  Doppler  shifts  and 
time  delays  is  correlated  with  a  reference  signal  containing  no  Doppler  shift  or  time  delay.  To 
understand  the  ambiguity  funetion  it  is  useful  to  start  by  introdueing  the  output  of  the  matehed 
filter  as: 

(3.1) 

In  Equation  3.1,  is  the  “reeeived  signal”,  s{t)  is  the  transmitted  or  “reference 
signal”,  s  *  (t)  is  the  eomplex  conjugate  of  the  reference  signal,  and  Tj^  is  the  estimated  time 
delay.  Using  complex  signal  representations,  the  transmitted  signal  ean  be  written  as 
I  u{t)  I  ,  where  u{t)  is  the  “complex  modulation  function”  whose  magnitude  is  the 

envelope  of  the  real  signal,  and  /g  is  the  carrier  frequency.  The  received  signal  is 
assumed  equal  to  the  transmitted  signal  exeept  for  a  Doppler  shift  of  and  a  time  delay 
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equaling  the  true  time  delay  T^.  Therefore,  the  matehed  filter  response  of  Equation  3.1 
beeomes: 

=  j  -  7;  exp[72;r(/o  +  f,  ){t  -  )]  exp[-72;^o  {t  -  )]dt  (3 .2) 

For  simplieity  in  understanding  this  equation,  the  origin  of  the  ambiguity  surfaee  is  taken 
to  be  the  true  time  delay  and  transmitted  frequeney,  i.e.,  =0,/^  =0  and  T^-Tj^  =  -T^  =  . 

After  some  simplifieations,  the  matehed  filter  output  beeomes: 

T 

ZiTj,  ,fd)  =  \  u{t)u  *  (t  +  )  QXT^{j2;tfj)dt  (3.3) 

0 

Here,  \Zi'rR,fd)\  will  be  ealled  the  ambiguity  function.  Properties  of  the  ambiguity 
funetion  inelude: 

1.  The  maximum  value  of  the  ambiguity  funetion  oeeurs  at  the  origin  as  expressed  by 
Equation  3.4.  This  oeeurs  sinee  the  reeeived  signal  arrives  at  the  reeeiver  with  the  same 
time  delay  and  zero  Doppler  shift  as  the  referenee  waveform;  this  is  effeetively  the  basie 
auto-eorrelation  response  of  the  waveform. 

Maximum  Value:  |  |Lx=l  =  i^Ef  (3.4) 
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2.  As  expressed  by  Equation  3.5,  the  ambiguity  surfaee,  ereated  by  plotting  the 
ambiguity  funetion  as  a  funetion  of  time  delay  and  Doppler  ,  is  symmetrie  along 
both  the  time  delay  and  frequeney  axes. 

Symmetry  Relation;  |  Zi-Tn-f,)  |"=|  ZiTnJd)  l"  (3.5) 

3.  The  total  volume  under  the  ambiguity  surfaee  is  eonstant  and  equals  (2E)  per 
Equation  3.6. 

Volume  Under  Surfaee:  jj|  =  {2EY  (3.6) 

The  maximum  value  property  of  Equation  3.4  states  that  the  value  at  the  origin  of  the 
ambiguity  funetion  is  {2EY  ,  where  E  is  the  total  Energy  eontained  in  the  signal.  Eikewise, 
Equation  3.6  indieates  that  the  total  volume  under  the  ambiguity  surfaee  is  also  {2EY  .  This  is  a 
total  volume  eonstraint  that  may  be  satisfied  by  many  different  surfaees  as  determined  by  the 
waveform;  whereas,  the  maximum  value  is  finite  at  the  origin. 

There  are  speeifie  methods  for  ealeulating  the  ambiguity  funetion.  The  first  teehnique 
eonsidered  multiplies  u(t)  with  u(t+Tj^  and  then  uses  an  Inverse  East  Eourier  Transform  (lEET) 
to  ealeulate  Equation  3.3.  The  eomplex  exponential  representation  of  the  signal,  as  used  in  this 
equation,  is  the  reason  why  lEET  teehniques  ean  be  used. 

In  the  seeond  technique,  for  a  given  frequency  Doppler  shift  ,  the  term 

u*  {t)Qx^{j27fj)  is  first  calculated.  Then,  u{t)  is  correlated  with  u*  {t)QXY){j27tfj) .  The 
correlator  output  for  different  values  yields  the  ambiguity  function  of  Equation  3.3. 
Correlation  in  the  time  domain  corresponds  to  multiplication  in  the  frequency  domain,  i.e.,  the 
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Fourier  spectrum  of  the  resultant  ambiguity  function  can  be  obtained  through  multiplication  of 
the  Fourier  transforms  of  u{t)  and  u*  .  Therefore,  one  specific  implementation 

of  the  second  technique  involves  calculating  the  Fourier  transforms  of  u{t)  and 
u*  {t)Qx^{j27fj)  for  different  values  first.  Then,  the  frequency  domain  transforms  are 

multiplied  and  the  result  is  the  Fourier  transform  of  the  ambiguity  function.  Hence,  by  taking  the 
IFFT  of  the  result  the  desired  response  of  Equation  3.3  is  obtained.  If  implemented  for  every 
Doppler  shift  /^,  the  ambiguity  surface  can  be  created  as  in  Figure  3.1.  This  method  is 
employed  for  this  research. 


0,8 


c 

•ri  0.6 


0 


Figure  3.1.  Normalized  Ambiguity  Surface  Plot  for  Rectangular  Sinusoidal  Pulse  Carrier 
Frequency  =  40A://zand  Pulse  Duration  =  0.0025  sec 
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Under  ideal  eonditions,  the  ambiguity  surfaee  eonsists  of  a  single  peak  having 
infinitesimal  thiekness  at  the  origin  and  zero  elsewhere,  i.e.,  an  impulse  funetion  having  no 
ambiguities  in  range  or  Doppler.  However,  beeause  of  ambiguity  funetion  properties  and  the 
resulting  surfaee  plot  it  is  not  possible  to  aehieve  the  ideal  impulse  funetion  at  the  origin.  The 
peak  width  along  the  time  delay  axis  determines  the  aeeuraey  and  resolution  of  time  domain 
measurements,  e.g.,  target  range;  and,  the  width  along  the  Doppler  axis  determines  the 
aehievable  aeeuraey  and  resolution  in  frequeney  domain  measurements,  e.g.,  veloeity.  In  both 
domains  a  smaller  peak  width  results  in  better  resolution  performanee.  Also,  the  higher  the  peak 
value  at  the  origin,  relative  to  the  surrounding  values,  the  better  the  measurement  aeeuraey 
beeomes. 

There  is  a  lot  to  learn  by  analyzing  the  ambiguity  surfaee  of  simple  waveforms  that  one 
ean  earry  into  analyzing  more  eomplex  waveforms.  By  examining  the  reetangular  sinusoidal 
pulse  ambiguity  surfaee  of  Figure  3.1  one  sees  a  triangle-shaped  response  of  the  ambiguity 
funetion  along  the  time  axis  where  there  is  no  Doppler  shift.  This  response  is  shown  in 
Equation  3.7,  whieh  is  the  evaluation  of  Equation  3.3  using  a  Doppler  frequeney  shift  of  zero. 

T 

=  +  (3.7) 

0 

The  expression  in  Equation  3.7  ean  also  be  seen  as  the  autoeorrelation  of  eomplex 
funetion  u{t)  .  The  time  axis  properties  of  Eigure  3.1,  along  with  the  ealeulations  from  Equation 
3.7,  yields  a  metrie  for  evaluating  the  time  domain  performanee  of  the  waveform,  the  simple 
reetangular  sinusoidal  pulse  in  this  ease.  In  the  eontext  of  this  researeh,  time  domain 
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performance  is  defined  as  the  achievable  measurement  accuracy  and  resolution  provided  by  the 
waveform. 
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Figure  3.2.  Normalized  Ambiguity  Surface  Plot,  Rectangular  Sinusoidal  Pulse  Carrier 
Frequency  =  40kHz  and  Pulse  Duration  =  0.005  sec  . 

Figure  3.2  is  the  ambiguity  surface  for  a  longer  sinusoidal  pulse  than  used  to  generate 
Figure  3.1,  t^=0.005sec  versus  =  0.0025sec .  The  sinusoid  amplitudes  are  equal.  The 

resultant  peak  at  the  origin  is  wider.  In  fact,  the  non-normalized  peak  response  at  the  origin  of 
Figure  3.2  is  higher  than  the  peak  response  of  the  shorter  pulse  ambiguity  surface  of  Figure  3.1, 
as  dictated  by  Equation  3.4.  Although  target  detection  capability  is  a  function  of  the  absolute 
peak  of  the  ambiguity  surface,  the  measurement  accuracy  and  resolution  are  not.  In  fact,  time 
delay  measurement  resolution  and  accuracy  is  a  function  of  frequency  signal  bandwidth  B. 
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Therefore,  for  time  domain  analysis,  one  ean  conclude  that  shorter  time  sinusoidal  pulses  provide 
better  range  measurement  accuracy  and  resolution. 

Ambiguity  surface  analysis  of  the  simple  sinusoidal  pulse  can  be  continued  by  analyzing 
Equation  3.3  along  the  frequency  axis  where  time  delay  is  zero.  The  resultant  equation  is  given 
in  Equation  3.8  below. 


ZiO,fd)  =  j w(Ow  *  it)Qxpij27tfj)dt  (3.8) 

0 

Equation  3.8  is  constructed  much  like  an  lEET  of  the  product  u(t)u  *(t)  .  Erom  Eourier 
transform  theory,  this  is  equivalent  to  correlation  of  the  lEET’s  of  u{t)  and  u*{t)  or  the 
autocorrelation  of  the  lEET  of  u{t) .  Erom  duality,  the  Eourier  transform  of  u{t) ,  being  a 
sinusoidal  pulse,  has  a  sine-shaped  lEET.  The  sine-shaped  lEET  of  u{t)  is  thus  correlated  with 
itself  yielding  the  zero  time  delay  surface  cut  of  the  ambiguity  surface. 
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Figure  3.3.  Zero  Time  Delay  Doppler  Frequeney  Axis,  Reetangular  Sinusoidal  Pulse  Carrier 
Frequeney  =  40^//z  and  pulse  duration  =  0.0025  see 


Figure  3.3  and  3.4  are  the  zero  time  delay  Doppler  frequeney  axis  plots  for  the  short  and 
long  sinusoidal  pulses,  respeetively. 

As  observed  by  eomparing  the  two  figures,  to  aehieve  good  resolution  in  Doppler 
frequeney  measurements  the  peak  width  at  the  origin  should  be  as  narrow  as  possible.  The  longer 
the  sinusoidal  pulse  duration  the  narrower  the  sine-squared  shape  beeomes  at  the  origin. 
Therefore,  the  Doppler  frequeney  measurement  aeeuraey  is  a  funetion  of  the  sinusoidal  pulse 
duration. 
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Figure  3.4.  Zero  Time  Delay  Doppler  Frequeney  Axis,  Reetangular  Sinusoidal  Pulse  Carrier 
Frequeney  =  40kHz  and  pulse  duration  =  0.005  see 

When  eonsidering  improved  measurement  performanee  in  both  the  time  and  Doppler,  a 
trade-off  relation  emerges  for  the  simple  reetangular  sinusoidal  pulse.  Namely,  to  improve  time 
domain  measurement  performanee  the  waveform  has  to  have  a  larger  bandwidth.  However,  to 
improve  Doppler  measurement  performanee  longer  pulses  are  required,  whieh  means  smaller 
signal  bandwidth  for  simple  sinusoidal  pulse.  This  trade-off  leads  to  introduetion  of  the  time- 
bandwidth  produet,  BT,  as  a  metrie  for  eharaeterizing  overall  waveform  performanee  in  a  radar 
or  PCL  system.  For  the  simple  baseband  reetangular  sinusoid  just  eonsidered  the  time-bandwidh 
produet  is  one. 
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3.3  Ambiguity  function  analysis  of  the  DAB  waveform 


The  analysis  of  more  eomplex  waveforms  is  not  as  straightforward  as  the  simple 
sinusoidal  pulse  ease.  In  ease  of  a  DAB  the  signal  eonsists  of  multiple  symbols  transmitted 
sequentially.  Each  symbol  consists  of  several  subpulses  having  pseudo-random  DQPSK 
modulation  that  are  eoherently  added  and  simultaneously  transmitted.  Details  of  the  DAB 
waveform  strueture  are  given  in  Chapter  2.  For  ambiguity  surface  analysis  of  the  DAB 
waveform  the  complex  u{t)  funetion  ean  be  formed  using  any  number  of  subpulses  and  DAB 
symbols.  The  total  number  of  DAB  symbols  used  in  forming  u{t)  determines  the  “eoherent 
integration  time”  for  ealeulations  in  Equation  3.3.  Coherent  integration  time  as  well  as  the 
number  of  subearrier  pulses  per  symbol  affects  DAB  waveform  performance  in  PCE  system 
applieations.  The  time-bandwidth  produet  of  a  DAB  PCE  waveform  is  generally  mueh  higher 
than  one.  In  the  DAB  waveform  case  the  resolution  of  time  delay  measurements  and  frequency 
shift  measurements  can  be  dealt  independently,  i.e.,  both  can  be  improved  simultaneously. 
Concerning  the  DAB  waveform  strueture  defined  as  COFDM,  simultaneous  improvement  is 
provided  by 

1.  Using  more  orthogonal-subearriers  per  symbol  to  effectively  inerease  the 
signal  bandwidth,  thereby  improving  time  delay  measurement  performanee, 
and 

2.  Increasing  the  coherent  integration  time  by  increasing  number  of  symbols  in 
the  representative  radar  waveform,  thereby  extending  the  signal  time  duration 
and  improving  Doppler  frequeney  measurement  performance  without 
degrading  time  domain  performance. 
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Note  that  increasing  the  number  of  symbols  per  coherent  integration  time  does  not  affect 
the  DAB  signal  bandwidth  -  DAB  signal  bandwidth  is  a  function  of  the  number  of  orthogonal 
subcarriers  per  symbol  and  the  time  duration  of  each  symbol. 

Two  methods  for  calculating  the  ambiguity  function  were  introduced  earlier.  For  this 
research,  the  second  method  has  been  adopted  and  implemented  per  the  following©(where  did 
you  introduce  it) 

1.  FFTs  of  u{t)  and  u  *  {t)Qx^{j27tfj)  for  the  desired  value  is  computed, 

2.  The  resultant  frequency  responses  are  multiplied  to  produce  the  frequency  domain 
representation  of  the  ambiguity  function  as  given  by  Equation  3.3, 

3.  An  IFFT  is  taken  to  yield  a  range  profile  (time  response)  at  the  selected  value. 

This  procedure  is  repeated  along  the  frequency  shift  axis  for  different  values  to  yield 
the  complete  ambiguity  function. 

For  any  given  signal  the  ambiguity  function  provides  a  measure  of  correlation  between 
reference  signal  M(t)and  the  Doppler-shifted  version  of  the  reference  signal  expressed  as 

u  *  .  In  fact,  the  reference  signal  itself  may  be  Doppler-shifted  as  well.  In  this 

case,  fj  axis  of  the  ambiguity  surface  represents  the  relative  Doppler  frequency  shift  between 
the  reference  u{t)  and  u  *  (t)exp(72;5f^t) ,  which  can  be  thought  as  the  return  of  the  reference 
signal  from  a  moving  object,  called  a  target. 

In  general,  any  waveform  used  for  calculating  an  ambiguity  function  perfectly  correlates 
with  itself  yielding  a  maximum  function  value  at  zero  time  delay  and  Doppler  frequency  shift. 
The  correlation  output  decreases  when  a  time  delayed  or  Doppler  shifted  version  of  the  reference 
is  correlated  with  the  reference  waveform.  This  relationship  can  be  used  for  determining 
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waveform  target  deteetion  eapability  through  ealeulation  of  the  ambiguity  funetion  and  plotting 
of  the  ambiguity  surfaee. 

In  PCL  applieations,  and  radar  applieations  in  general,  the  time-delay  and  Doppler  shift 
of  a  reeeived  target  return  is  unknown.  Assuming  an  independent  referenee  signal  ean  be  ereated 
at  the  reeeiver  one  simple  parametrie  deteetion  and  estimation  teehnique  for  unknown  target 
eharaeteristies  eonsists  of: 

1 .  Adding  a  time  delay  and  Doppler  shift  to  the  referenee  waveform, 

2.  Correlating  this  waveform  with  the  reeeived  target  return  signal, 

3.  Repeating  steps  1  and  2  until  the  entire  “deteetion  surfaee”  is  generated, 

4.  Loeating  the  peak  eorrelation  response  on  the  surfaee,  and 

5.  Estimating  the  unknown  target  parameters  as  those  eoineiding  with  the  time  delay 
(range  estimate)  and  Doppler  (relative  veloeity  estimate)  values  where  the  peak 
oeeurs. 

Theoretieally,  the  highest  eorrelation  response  oeeurs  when  both  eomponents  entering  the 
eorrelation  are  equal.  Another  possible  outeome  is  that  a  peak  value  is  not  suffieient  to  trigger 
target  deteetion,  e.g.,  a  target  response  is  present  but  it  has  insuffieient  energy  to  produee  a 
surfaee  response  above  the  noise  floor.  Alternately,  if  a  target  is  not  present,  the  ealeulated 
ambiguity  funetion  will  always  have  a  peak  response  somewhere.  Whether  or  not  a  peak 
response  eorresponds  to  a  target  or  not  is  ultimately  determined  by  the  signal  proeessing 
teehniques  used.  Ambiguity  funetion  ealeulation  and  ambiguity  surfaee  analysis  of  the  referenee 
waveform  response  is  a  vital  input  for  any  signal  proeessing  effort  for  deteeting  targets  with  a 
matehed  filtering  approaeh.  For  example,  thresholding  ean  utilize  a  peak  to  sidelobe  ratio  of 
ambiguity  surfaee. 
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The  real  utility  of  the  ambiguity  surfaee  ean  be  best  exploited  if  the  inputs  to  the 
ambiguity  funetion  formula  of  Equation  3.3  are  realistie,  i.e.,  when  the  referenee  waveform  is 
being  Doppler-shifted,  the  signal  return  of  the  referenee  waveform  from  a  target  with  a  relative 
veloeity  eorresponds  to  the  applied  Doppler-shift.  This  implementation  is  relatively  easy  for 
simple  sinusoidal  pulsed  waveforms.  If  the  referenee  waveform  eonsists  of  more  than  one 
sinusoidal  subpulse,  as  the  ease  in  DAB,  the  situation  gets  mueh  more  eomplex.  Namely,  eaeh 
subpulse,  or  DAB  subearrier,  experienees  a  different  Doppler  frequeney  shift  as  dietated  by  its 
earner  frequeney  (or  wavelength)  and  the  relative  veloeity  of  a  moving  target.  DAB  Mode  I  has 
1536  different  subearrier  eenter  frequeneies  [14].  Therefore,  shifting  all  subearriers  by  the  same, 
eonstant  Doppler  frequeney  is  not  a  realistie  approaeh.  However,  there  is  an  advantage  gained 
by  applying  a  eonstant  frequeney  shift  to  all  subearriers.  Namely,  the  ambiguity  surfaee 
ealeulations  beeome  mueh  faster  sinee  a  simple  IFFT  ean  be  used  to  ereate  the  Doppler  shifted 
signal. 

This  next  seetion  explains  how  a  simple  IFFT  ean  be  used  to  ereate  the  frequeney  shifted 
DAB  waveform  (assuming  a  eonstant  Doppler  shift  for  all  subearriers)  as  well  as  an  explanation 
of  the  theoretieally  eorreet,  but  more  tedious  way  to  ealeulate  the  frequeney  dependent  Doppler 
shifted  DAB  signal. 

The  DAB  Mode  I  waveform  eonsists  of  1536  orthogonal  subearriers  of  1.246 
milliseeonds  duration  transmitted  simultaneously  as  single  DAB  symbol.  As  mentioned  in 
Chapter  2,  there  is  a  eyelie  prefix  added  to  the  beginning  of  eaeh  subearrier,  i.e.,  0.246 
milliseeonds  from  the  end  of  eaeh  subearrier  is  copied  and  pasted  at  the  beginning  [14]. 
Basically,  each  DAB  symbol  can  be  written  as  the  sum  of  1536  orthogonal  sinusoids  as  given  in 
Equation  3.9.  Detailed  information  related  to  DAB  waveform  is  included  in  Chapter  2. 
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(3.9) 


1536 

1=1 

Given  a  DAB  signal  is  transmitted,  the  return  signal  from  a  moving  target  will  induee  a 
different  Doppler  frequency  shift  on  each  subcarrier  as  determined  by  the  subpulse  center 

frequency  /  .  The  target  return  will  be  similar  to  the  expression  in  Equation  3.10; 

1536 

^  cos[(/;.  +  )lm  +  6>,  ]  (3.1 0) 

1=1 

Equation  3.10  is  the  proper  analytic  form  of  a  Doppler  shifted  DAB  symbol,  clearly 
indicating  a  subpulse  carrier  dependent  Doppler  component.  To  accurately  simulate  this  Doppler 
shifted  symbol  one  has  to  induce  a  different  Doppler  shift  on  each  subcarrier  frequency  for  a 
given  relative  velocity.  Eor  computer  simulations  this  approach  is  impractical.  It  is  desirable  to 
have  a  fast  algorithm  for  computing  the  Doppler  shifted  DAB  signal  without  losing  important 
information.  One  fast  approach  is  to  assume  a  constant  Doppler  shift  (subcarrier  frequency 
independent)  for  every  subcarrier,  thereby,  that  fast  lEET  technique  can  be  used  to  create  target 
return  signals.  To  demonstrate  this  approach.  Equation  3.10  is  rearranged  and  rewritten  to 
include  a  constant  Doppler  shift  as  illustrated  in  Equation  3.1 1.  By  applying  some 

trigonometric  identities.  Equation  3.12  through  Equation  3.14  show  the  fast  lEET  creation  of  a 
Doppler  shifted  DAB  signal. 

1536 

Sj  =  Y,cos[{f,2m  +  e,)  +  f,2m]  (3.11) 

/=! 

1536 

St  =  '^  cos(/^  2m)  cos(/.  2m  +  9.)-  sin(/^  2m)  sin(/.  2m  +  9.)  (3.1 2) 

1=1 
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(3.13) 


1536  1536 

Sj.  =  cosif, 2m)Y,cosif  2m  +  9,)- sin(/^ 2m)Y,  sin(/;. 2m  +  6>, ) 

i=\  i-1 

Sj  =cos{f^2m)KQ{IFFT{g{e,))\-sm{f^2m)\m{IFFT{g{e.))\  (3.14) 

1536 

In  going  from  Equation  3.13  to  Equation  3.14,  the  ^  cos{f.2m  +  0.)  term  is  replaced 

/=1 

1536 

by  RQ[IFFT(g(9J)]  and  the  ^  sin(/.  )  term  is  replaced  by  lm[IFFT(g(9J)] .  This 

1=1 

lEET  representation  was  introduced  in  Chapter  2  to  facilitate  the  efficient  creation  of  DAB 
subcarriers.  Therefore,  under  the  assumption  of  constant  Doppler  shift  for  all  subcarriers,  the 
lEET  can  be  used  for  DAB  signal  formation  and  for  constructing  a  Doppler  shifted  target  return. 
Clearly,  assuming  a  constant  subcarrier  Doppler  shift  produces  error.  The  question  is,  how  much 
error  is  introduced.  The  amount  of  error  created  by  this  assumption,  and  its  effect  on  ambiguity 
function  calculation  and  surface  generation,  is  demonstrated  in  Chapter  4. 

3.4  Practical  issues  with  DAB-PCE  system 

This  section  introduces  a  detrimental  factor  present  in  the  implementation  of  a  realistic 
PCE  system,  namely,  the  presence  of  a  strong  reference  signal  component  in  the  target  channel 
of  a  PCE  receiver.  A  PCE  receiver  dedicated  to  exploiting  DAB  waveform  may  use  matched 
filter  techniques  for  signal  processing  and  target  detection.  Eor  this  particular  processing 
technique,  at  least  two  receiver  channels  are  necessary.  One  channel  is  dedicated  to  receiving 
the  reference  waveform,  defined  as  the  direct  response  from  a  line-of-sight  (EOS)  transmitter. 
The  second  receiver  channel  called  the  target  channel,  i.e.,  the  antenna  of  this  receiver  channel  is 
specifically  directed  towards  a  surveillance  region  to  receive  target  returns  (echoes)  of  the 
reference  waveform.  It  is  practically  impossible  to  completely  isolate  the  two  receiver  channels. 
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i.e.,  one  component  of  the  reference  channel  signal  is  a  weaker  target  return  and  one  component 
of  the  target  channel  signal  is  a  stronger  reference  signal  response. 

The  target  return  component  in  the  reference  channel  does  not  pose  a  serious  problem 
since  this  component  is  generally  very  much  lower  in  power.  Thus,  it  does  not  affect  system 
performance.  However,  reference  signal  leakage  into  the  target  channel  signal  significantly 
affects  system  performance.  A  rigorous  mathematical  analysis  of  this  detrimental  effect  is 
provided  in  reference  1 1  with  important  conclusions  from  this  analysis  provided  below. 

Equation  3.15  and  Equation  3.16  are  the  received  signals  for  both  channels,  including  the 
interfering  components. 

=  +  (3.15) 

w, (0  =  Pau{t)  +  /3^u{t  -  h ) (3.16) 

In  Equation  3.15  u^{t)  is  the  received  signal  in  the  reference  channel  and  contains  two 
components,  reference  signal  component  a^uit)  and  target  return  component 
a^uit  -  ti)exp{j27tfit)  .  The  target  return  component  is  nothing  more  than  the  reference  signal 
delayed  in  time  by  h  (a  function  of  target  range)  and  shifted  in  frequency  by  /j  (a  function  of 
target  relative  velocity).  As  shown  in  Equation  3.16,  the  received  signal  in  the  target  channel 
u,{t)  contains  similar  components  with  different  weights  of  and  If  the  weights  are 

compared,  one  sees  that  the  largest  weight  is  ,  the  weight  of  the  reference  component  in  the 

reference  channel.  Although  the  target  channel  is  dedicated  to  receiving  and  processing  the 
target  return  signal  and,  therefore,  it  is  desirable  for  the  target  return  component  to  dominate  in 
this  channel;  in  practice  there  is  often  strong  leakage  from  the  EOS  transmitter  and  the  reference 
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component  having  weight  dominates  in  the  target  channel.  The  weakest  link  is  the  target 
return  component  in  the  reference  channel  having  weight  (assuming  the  target  channel  gain  is 
higher  than  the  reference  channel  gain  in  the  target  direction). 

When  both  receiver  channel  signals  are  correlated  with  each  other,  the  output  will  be 
given  as  expressed  by  Equation  3.17: 

^  MF  ~  ^Zi^R  ^Zi^R  ~h^fd  ~  /i  )  “*“  ^Zi^d  '^hifd'^fl)  (3-17) 

a  =  +  a.fil  QXY>{j27f,nty  -  j27tf,T^) 

b  =  ao/3*  Qxp{-j2;tf^Tj,)  (3.18) 

c  =  Qxpijti  {27tf^n  +  27tf,)) 

If  Equation  3.17  is  plotted  as  a  correlation  surface  the  three  peak  responses  described  by 
Equation  3.18  become  visible.  The  first  peak  response  occurs  at  the  origin  and  is  the  first  term  in 
Equation  3.17  having  weight  a.  In  this  case,  the  reference  signal  components  in  the  target  and 
reference  channels  correlate  well  and  result  in  a  peak  response  at  the  origin  having  weight  . 
At  the  same  time,  the  target  return  components  in  eaeh  channel  correlate  well  with  each  other 
resulting  in  a  peak  response  having  weight  .  The  second  peak  is  located  at  {-ti ,  -/i)  and 
has  weight  CC^Pl.  This  ‘spurious’  peak  results  from  the  cross-correlation  of  the  target  return  in 
reference  channel  with  reference  signal  leakage  in  the  target  channel.  A  third  peak  is  present  at 
(tj  ,  /j)  and  is  appropriately  located  with  weight  .  This  peak  is  the  cross-correlation 

response  of  the  reference  component  in  the  reference  channel  with  the  target  return  component  in 
the  target  channel. 
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The  detection  challenge  lies  in  distinguishing  the  correct  (desired)  peak  response  from  the 
other  two.  The  peak  located  at  -ti  and  - /j  has  the  smallest  weight  and  is  generally  not 

harmful.  However,  the  peak  at  the  origin  is  weighted  much  more  heavily  )  and  provides  a 

significant  dominating  response.  To  be  successful  in  target  detection,  any  signal  processing 
architecture  must  be  designed  to  eliminate  this  peak. 

In  reality  other  components  not  previously  mentioned  may  be  present  in  the  receiver 
target  channel,  namely,  the  multipath  response(s)  of  the  reference  signal.  Details  concerning  the 
multipath  response  in  DAB  SFN  architectures  are  covered  in  Chapter  2.  The  multipath 
component(s)  can  significantly  impact  PCL  system  performance.  Specifically,  this  component 
may  be  stronger  than  the  actual  target  return,  causing  additional  peaks  in  the  correlation  surface 
as  plotted  using  Equation  3.17.  Therefore,  the  direct  signal,  defined  as  the  reference  and 
multipath  components  in  receiver  target  channel,  should  be  removed.  Ideally,  this  enables  the 
target  return  to  trigger  a  visible  peak  in  the  correlation  surface,  resulting  in  reliable  target 
detection  and  parametric  measurement  (estimation)  of  range  and  relative  velocity. 

This  research  proposes  an  algorithm  for  eliminating  the  direct  signal,  although  there  are 
no  claims  of  optimality  or  promise  of  complete  removal.  The  basic  idea  is  to  first  remove  or 
attenuate  the  direct  signal  component  in  the  target  channel  prior  to  detection  using  a  matched 
filtering  process.  The  basic  idea  behind  the  technique  developed  for  this  research  involves  signal 
projection  via  a  linear  algebra  matrix  operation,  i.e.,  projection  of  vectors  into  column  space  of 
matrices  [22,23],  as  detailed  in  Equations  3.19  through  3.21. 

P  =  A{A^Ay^A^  (3.19) 
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R  =  I-P 


(3.20) 


v  =  Rs 


(3.21) 


In  Equation  3.19,  matrix  P  projects  a  given  vector  into  the  column  space  of  A,  where 
columns  of  A  may  be  selected  as  orthogonal  components  of  the  direct  signal  component  in  the 
target  channel.  Consequently,  P  projects  any  given  vector  (samples  of  the  received  signal)  into 
the  space  defined  by  orthogonal  components  of  the  direct  signal. 

In  Equation  3.20,  identity  matrix  /  has  the  same  dimension  as  P.  Vector  R  of 
Equation  3.21,  if  post-multiplied  with  a  vector  (samples  of  received  signal  in  target  channel), 
projects  the  vector  into  the  null  space  of  matrix  A  columns.  If  columns  of  ^4  represent  orthogonal 
components  spanning  the  direct  signal  space,  assuming  s  in  Equation  3.21  is  a  column  vector  of 
received  signal  samples  in  the  target  channel,  v  will  be  the  projection  of  s  into  a  space  orthogonal 
to  the  direct  signal.  In  ideal  cases,  i.e.,  A  can  be  calculated  error-free,  this  operation  yields  signal 
V,  which  is  completely  void  of  direct  path  components. 

A  =  [sls^^sl....s:]  (3.22) 

In  equation  3.22,  is  the  vector  containing  the  time  samples  of  subcarrier  x  waveform 
and  n  is  the  number  of  subcarriers  per  DAB  symbol.  The  received  signal  s  in  Equation  3.21  can 
be  written  as  the  sum  of  direct  signal  target  return  noise  component  ^^^i^^as  in 

Equation  3.23. 
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^direct  ^target  ^ noise 


(3.23) 


^  direct 


can  be  written  as  the  sum  of  the  subcarrier  time  samples  (Equation  3.24) 


1,2,3,  it, 
S  direct  =Sc+S,  +S^  +....  +  5, 


(3.24) 


Finally  Equation  3.21  becomes  as  the  following: 


v  =  (l-p).  ) 


(3.25) 


v  =  (I-P).(sl+  5"  +  +....  +  5;  )  +  (/  -  E)  •  ) 


(3.26) 


The  first  term  of  Equation  3.26  is  zero  because  (/  -  P)  spans  a  space  orthogonal  to  , 
comprised  of  the  sum  of  reference  signal  and  multipath  components  as  explained  in  Chapter  2. 
Therefore,  resultant  v  void  of  . 

Two  issues  must  be  addressed  with  regard  to  this  scheme: 

1.  The  orthogonal  components  spanning  the  direct  signal  space  cannot  be 
calculated  error-free. 

2.  The  target  return  component  in  received  signal  s  will  experience  attenuation 
during  the  projection  process. 

These  two  issues  are  discussed  and  analyzed  in  greater  detail  in  Chapter  4. 

3.5  Implementation  of  direct-path  filter 

In  implementing  the  algorithm  and  analyzing  its  performance  the  first  step  taken  included 
creating  the  received  signals.  The  received  signal  in  both  channels  is  assumed  to  have  four 
components-  a  reference  signal,  a  multipath  signal  (together  with  the  reference  signal  these  two 
are  called  the  direct  signal),  a  target  return  and  channel  noise. 
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The  reference  signal  is  created  first  in  the  process  followed  by  target  return  signal 
generation  using  a  user-defined  relative  Doppler-velocity.  The  input  relative  velocity  is  used  for 
inducing  correct  Doppler  shift  on  individual  DAB  subcarriers.  Next,  the  multipath  signal  is 
created  using  a  user-defined  path  delay  relative  to  the  reference  signal.  The  target  return  signal 
power  and  multipath  signal  power  are  adjusted  according  to  user-defined  values  expressed  in  dB. 
The  channel  noise  component  is  assumed  Gaussian  and  it  power  is  adjusted  to  provide  the  proper 
SNR  value. 


Figure  3.5.  Block  Diagram  of  Direct  Path  Filter  signal  processing 
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As  explained  in  Chapter  2  the  DAB  signal  is  transmitted  in  frames.  The  first  non-zero, 
data  modulated  symbol  of  every  DAB  frame  transmission  is  the  phase  referenee  symbol.  A 
phase  referenee  symbol  having  a  predetermined  eode  sequenee  is  known  a  priori  by  the  reeeiver 
and  used  to  determine  subearrier  phase  thresholds  (this  is  required  for  DQPSK  demodulation). 
This  represents  an  important  part  of  the  DAB  waveform  demodulation  proeess,  where  indireetly, 
the  channel  frequency  response  is  estimated  and  used  to  achieve  “near  error-free”  data 
demodulation. 

Two  different  phase  reconstruction  cases  have  been  simulated  for  this  research.  The  first 
case  assumes  perfect  construction  of  the  phase  thresholds  for  each  subcarrier,  which  requires  the 
phase  reference  symbol  to  be  made  up  of  only  the  direct  signal  component.  In  this  case,  the 
received  signal  can  be  demodulated  correctly.  Although  in  the  first  case,  the  received  signal  for 
each  symbol  (except  the  phase  reference  symbol)  has  a  target  return  and  noise  component  in 
addition  to  direct  path,  demodulation  based  on  the  perfect  phase  reference  symbol  and  then 
remodulating  the  resulting  data  with  the  phase  values  found  from  the  perfect  phase  reference 
symbol  creates  a  perfect  base.  The  perfect  base  can  be  defined  as  the  direct  path  only  component 
of  the  received  signal  calculated  by  demodulating  the  received  signal  through  the  perfect  phase 
reference  signal  and  then  remodulated  again.  The  subcarriers  that  make  the  perfect  base  are  used 
by  the  filter  algorithm  to  form  matrix  A  of  Equation  3.22. 

For  case  two,  more  realistic  conditions  are  considered.  Namely,  all  symbols,  including 
the  phase  reference  symbol,  suffer  degradation  due  to  noise.  Also,  the  target  return  component 
in  the  reference  channel  has  a  negative  effect  on  DAB  demodulation.  The  direct  signal 
component  in  conjunction  with  the  channel  noise  changes  the  demodulation  phase  thresholds  for 
each  subcarrier.  Consequently,  even  if  the  data  can  be  demodulated  correctly,  in  the  subsequent 
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re-modulation  process  for  creating  the  perfect  base  for  the  direct  path  filter,  the  calculated  phase 
values  will  not  accurately  represent  the  direct  path  component  for  each  subcarrier.  This  is  the 
result  of  a  noisy  phase  reference  symbol  causing  incorrect  phase  thresholds  for  each  subcarrier. 

In  the  proposed  direct  path  fdter  design,  the  phase  reference  symbol  is  demodulated  first. 
Phase  threshold  values  for  each  subcarrier  used  in  demodulation  of  DQPSK  waveform  are 
determined  as  a  result  of  demodulating  the  phase  reference  symbol. 

The  next  step  in  the  algorithm  uses  the  calculated  thresholds  from  the  phase  reference 
symbol  to  demodulate  the  data.  In  demodulation  the  phase  of  each  subcarrier  in  a  symbol  is 
compared  with  the  phase  thresholds.  After  demodulation,  the  data  is  re-modulated  to  the  correct 
phase  values  for  each  subcarrier  to  create  a  base  for  the  filtering  design. 

To  form  filtering  matrix  R,  as  defined  in  Equation  3.20,  each  subcarrier  waveform  after 
re-modulation  is  used  to  form  matrix  R,  which  is  assumed  to  span  the  direct  signal  space.  To  do 
this  time  samples  of  each  subcarrier  are  placed  in  a  column  of  the  matrix.  Given  the  subcarriers 
are  orthogonal  to  each  other,  the  columns  of  the  matrix  are  orthogonal  as  well.  The  resultant 
matrix,  denoted  A  in  Equation  3.19,  is  used  to  calculate  the  P  and  R  matrices  of  Equations  3.19 
and  3.20,  respectively.  Next,  R  is  multiplied  with  the  received  signal  per  Equation  3.21.  The 
resultant  signal  v  is  the  filter  output.  In  the  perfect  case  where  a  noise-free  phase  reference 
symbol  yields  a  perfect  base  the  signal  at  the  filter  output  is  void  of  any  direct  signal.  However, 
such  ideal  conditions  are  impractical. 

To  illustrate  the  noise  effects  and  characterize  filter  performance  of  the  proposed  design, 
detection  surfaces  are  created  for  the  received  signal  before  (pre-filtered)  and  after  filtering 
(post-filtered).  The  detection  surface,  similar  to  an  ambiguity  surface,  represents  the  correlation 
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of  a  reference  signal  with  the  received  signal  at  different  time  delay  and  Doppler  shift  values. 
Filter  performance  results  are  presented  in  Chapter  4. 

3.6  Summary 

This  chapter  introduced  the  ambiguity  function  as  a  tool  for  evaluating  waveforms  of 
opportunity  for  PCL  applications.  Ambiguity  function  properties  are  stated.  Next,  DAB  related 
ambiguity  function  issues  are  discussed  followed  by  discussions  on  the  effects  of  receiving  a 
direct  path  signal  in  the  target  channel.  The  performance  of  detection  using  correlation 
techniques  is  presented  as  well.  Finally,  a  unique  filtering  technique  is  proposed  to  eliminate  the 
detrimental  effects  of  direct  path  leakage  into  the  target  channel  by  removing  or  significantly 
attenuating  it  before  matched  filter  detection. 

The  next  chapter  includes  the  ambiguity  function  implementation  for  DAB  analysis  as 
implemented  in  MATLAB®.  It  also  provides  performance  results  for  the  proposed  filter  design 
specifically  created  to  eliminate  the  detrimental  effects  of  direct  path  leakage  using  a  MATLAB® 
implementation  of  the  filter. 
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CHAPTER  4 


SIMULATION  RESULTS  AND  ANALYSIS 


4.1  Introduction 

In  this  chapter  the  concepts  of  previous  chapters  are  used  to  conduct  performance 
analysis  of  a  Digital  Audio  Broadcast  (DAB)  signal  and  to  implement  a  filter  designed  to 
eliminate  the  reference  and  multipath  component  (together  called  direct  path)  of  the  received 
signal  in  the  target  channel  of  the  PCL  receiver.  This  chapter  starts  with  a  discussion  about  DAB 
related  issues  in  ambiguity  function  calculations.  Lollowing  that  there  is  an  evaluation  of  the 
performance  of  DAB  as  a  waveform  of  opportunity  for  PCL  systems.  After  the  performance 
evaluation,  a  filter  design  based  on  the  projection  operator  introduced  in  Chapter  3  is  created  and 
then  implemented  on  the  simulated  DAB-PCL  received  signal.  The  chapter  concludes  with  an 
analysis  of  the  filter  design  performance. 

4.2  Ambiguity  function  analysis  of  the  DAB  waveform 

The  discussions  in  Chapter  2  included  creation  of  a  DAB  waveform  through  the 
implementation  of  an  inverse  LET  algorithm.  This  algorithm  makes  DAB  waveform  generation 
very  fast  and  efficient.  As  mentioned  in  Chapter  3,  ambiguity  function  calculation  includes  the 
ereation  of  a  Doppler  shifted  signal.  Next,  a  Last  Lourier  Transform  (LET)  of  both  the  signal  and 
a  Doppler  shifted  version  of  the  signal  are  calculated  and  multiplied  together.  An  inverse  LET  of 
the  resultant  product  yields  the  waveform  ambiguity  response  for  the  caleulated  Doppler  shift  as 
a  function  of  time.  Repeating  this  step  for  each  desired  Doppler  response  is  necessary  to  form 
the  complete  waveform  ambiguity  surface  for  a  given  set  of  Doppler  frequencies. 
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One  major  challenge  in  forming  the  ambiguity  surface  for  a  DAB  waveform  via 
computer  modeling  and  simulation  exists  because  of  the  complex  structure  of  the  waveform.  A 
DAB  symbol  consists  of  1536  orthogonal  sinusoids  that  are  transmitted  simultaneously  for  a 
symbol  duration  of  1.246  ms.  Considering  the  induced  Doppler  frequency  of  a  moving  target  is 
calculated  according  to  Equation  4. 1 . 

fd=2*v*f,/c  (4.2) 

To  accurately  model  a  Doppler  shifted  DAB  signal,  one  has  to  calculate  the  appropriate 
Doppler  shift  for  every  subcarrier  individually.  Subcarrier  center  frequency  is  denoted  by  C, 
where  v  is  the  relative  speed  of  the  target  inducing  an  echo  signal  at  the  receiver  target  channel. 
The  individual  DAB  Doppler  subcarrier  shifts  are  not  the  same  for  a  given  relative  velocity 
because  of  the  difference  in  subcarrier  center  frequency.  Another  source  of  complexity  in 
creating  a  reliable  Doppler  shifted  DAB  signal  arises  from  the  bandwidth  of  individual 
subcarriers.  Signal  bandwidth  defines  a  range  of  frequencies  (frequency  bin)  for  each  sinusoid 
as  shown  in  Figure  4.1.  To  be  theoretically  correct,  one  has  to  calculate  induced  Doppler  for 
each  frequency  component  in  the  signal.  Practical  concerns,  together  with  the  fact  that  each 
sinusoid  has  a  narrowband  spectrum,  dictates  those  Doppler  calculations  to  be  made  only  for 
center  frequencies  of  each  sinusoid. 
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The  frequency  bin  of  the  subcarrier  with  f^.  as  center  frequency 


Figure  4.1.  DAB  Symbol  Bandwidth 


Figure  4.1  shows  basic  DAB  signal  structure  in  the  frequency  domain.  Sinusoids  that 
have  been  phase-shift  keyed,  called  subcarriers,  have  a  sine-shaped  frequency  spectrum.  Every 
subcarrier  has  distinct  center  frequency  that  defines  the  subcarrier.  All  subcarriers  have  an 
infinite  bandwidth  as  determined  by  the  sine-shape.  As  for  the  Doppler  shifted  DAB  signal 
calculations,  for  each  subcarrier  every  frequency  component  of  the  signal  is  affected  differently 
according  to  the  formula  of  Equation  1.  In  Eigure  4.1,  a  subcarrier  frequency  bin  is  defined  as 
the  signal  bandwidth  containing  90%  of  the  total  subcarrier  power.  Practically,  subcarrier 
Doppler  shift  is  calculated  based  on  the  center  frequency  of  each  bin.  Therefore,  the  signal’s 
spectral  shape  remains  unchanged  after  shifting.  This  fact  makes  analysis  of  Doppler  shifts 
easier.  Therefore,  assuming  there  are  1536  subcarriers  starting  with  a  center  frequency  of  1  kHz 
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at  baseband,  Doppler  frequeney  may  be  calculated  based  on  one  single  frequency  located  at  the 
center  of  the  total  DAB  signal  bandwidth.  This  approach  has  an  inherent  advantage  in  ambiguity 
function  calculation  of  the  DAB  waveform.  Specifically,  this  approach  enables  us  to  use  fast 
DAB  symbol  creation  algorithms  based  on  IFFT  and  Doppler  shifted  signals  are  easily  created 
(Chapter  3). 

The  issue  that  one  has  to  pay  attention  is  the  introduction  of  an  error  with  the  assumption 
of  constant  Doppler  frequency  for  all  of  the  subcarriers.  The  next  section  deals  with  visualizing 
the  error  caused  by  implementing  the  fast  DAB  waveform  generation  techniques  explained 
above,  techniques  which  have  been  readily  adopted  by  other  researchers  with  minimal 
justification  [18]. 

To  see  the  extent  of  the  problem,  two  ambiguity  surface  generation  techniques  are 
considered.  For  method  one,  a  constant  Doppler  frequency  shift  resulting  from  a  moving  target  is 
assumed  induced  on  each  subcarrier,  i.e.,  the  amount  of  Doppler  shift  induced  on  each  subcarrier 
is  independent  of  the  actual  subcarrier  frequency.  For  method  two,  the  ambiguity  surface  plot 
represents  the  correct  Doppler  shifted  signal  in  the  simulation,  i.e.,  in  creating  the  Doppler 
shifted  signal  the  induced  Doppler  resulting  from  a  moving  target  is  calculated  independently  for 
each  subcarrier,  the  subcarriers  are  then  added  to  form  the  composite  Doppler  shifted  signal. 

In  analyzing  three-dimensional  ambiguity  surface  plots  from  both  techniques,  it  is 
impossible  to  visually  discern  any  difference  between  the  approaches.  So,  to  do  a  comparison  of 
results  from  both  approaches,  the  response  along  a  principal  axis  of  the  surface  is  analyzed. 
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Figure  4.2.  Method  1:  Independent  Subearrier  Doppler,  Zero  Doppler  DAB  Ambiguity  Surface 
Cut  (30  Subcarriers  and  20  Symbol  Coherent  Integration  Time) 


Figure  4.2  shows  the  zero  Doppler  ambiguity  surface  cut  (time-delay  response)  for  the 
downscaled  DAB  waveform,  where  one  DAB  symbol  consists  of  30  subcarriers  and  the  coherent 
integration  time  consists  of  20  DAB  symbols.  The  response  is  normalized  and  expressed  in  dB. 
As  shown  in  Chapter  3,  the  response  also  represents  the  autocorrelation  of  the  signal.  In 
practice,  each  DAB  symbol  consists  of  1536  orthogonal  sinusoidal  subcarriers  (for  DAB  Mode 
I).  The  time  axis  extends  to  about  1.0  ms,  which  corresponds  to  a  monostatic  range  of  150  km. 
Range  sidelobes  in  this  case  appear  7  to  15  dB  below  the  maximum  response.  To  create  this 
figure,  the  Doppler  shifted  signal  is  calculated  using  the  inverse  FFT  technique  described  in 
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Chapter  3,  i.e.,  a  constant  Doppler  frequency  shift  is  applied  to  each  subcarrier  center  frequency, 
independent  of  individual  subcarrier  frequencies. 

Figure  4.3  is  the  same  reference  signal  response  except  that  method  two  has  been  used  to 
generate  the  response.  For  generation  using  method  two,  the  Doppler  shifted  signal  is  calculated 
to  form  the  ambiguity  function  as  follows:  each  subcarrier  has  a  different  amount  of  induced 
Doppler  due  to  having  different  center  frequencies  (Doppler  dependence  on  wavelength  is  thus 
incorporate).  This  approach  is  more  realistic  since  each  subcarrier  Doppler  shift  is  a  function  of 
the  frequency  bin  content  defined  by  the  subcarrier  center  frequency. 
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Figure  4.3.  Method  2:  Dependent  Subcarrier  Doppler,  Zero  Doppler  DAB  Ambiguity  Surface 
Cut  (30  Subcarriers  and  20  Symbol  Coherent  Integration  Time) 
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Visual  comparison  of  Figure  4.2  with  Figure  4.3  reveals  minimal  difference  in  the 
responses.  Therefore,  to  help  quantify  performance  of  the  two  ambiguity  surface  generation 
techniques,  Figure  4.4  is  presented  and  represents  the  ‘difference’  of  the  outputs  taken  at  each 
time  delay.  These  results  indicate  an  absolute  difference  of  less  than  lO'^"^  for  all  time  delays. 
Clearly,  both  approaches  yield  nearly  identical  outputs  for  the  case  being  considered. 
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Figure  4.4.  Difference  Plot  of  Zero  Doppler  Ambiguity  Surfaces  Generated  Using  Independent 

and  Dependent  Subcarrier  Induced  Doppler 


Next,  the  same  type  of  analysis  is  conducted  for  the  zero  time-delay  ambiguity  surface 
cut  (Doppler  response).  Figure  4.5  presents  results  for  method  one  where  Doppler  shifted  signal 
is  calculated  fast  using  IFFT. 
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Figure  4.5.  Method  1:  Independent  Subearrier  Doppler,  Zero  Time-Delay  DAB  Ambiguity 
Surface  Cut  (30  Subcarriers  and  20  Symbol  Coherent  Integration  Time) 


In  Figure  4.5,  the  induced  Doppler  shift  is  calculated  as  a  function  of  relative  target 
velocity  given  along  the  x-axis.  For  this  case,  every  subcarrier  is  affected  by  a  constant  Doppler 
frequency  shift  (independent  of  subcarrier  center  frequency)  as  a  function  of  a  chosen  center 
frequency,  e.g.,  the  chosen  frequency  can  be  the  carrier  center  frequency  in  the  middle  of  the 
DAB  signal  spectrum.  Therefore,  for  this  example  where  every  DAB  symbol  has  30  subcarriers, 
the  center  frequency  of  the  15*’^  subcarrier  is  chosen  as  the  frequency  to  base  Doppler  shift 
calculations  on  for  the  entire  symbol. 
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Figure  4.6.  Method  2;  Dependent  Subearrier  Doppler,  Zero  Time-Delay  DAB  Ambiguity 
Surface  Cut  (30  Subcarriers  and  20  Symbol  Coherent  Integration  Time) 


Figure  4.6  is  a  plot  of  the  zero  time-delay  ambiguity  surface  cut  generated  using  method 
two  (subcarrier  dependent  Doppler).  In  this  technique,  the  target  relative  velocity  is  used  with 
each  subcarrier  center  frequency  to  calculate  the  unique  Doppler  shift  for  that  subcarrier,  i.e.. 
Equation  4.1  is  used  with  each  subcarrier  center  frequency  to  determine  the  corresponding 
Doppler  frequency  shift.  In  both  Figure  4.5  and  Figure  4.6,  the  velocity  sidelobes  tend  to  be  10 
to  15  dB  below  the  maximum  response.  As  in  the  previous  zero  Doppler  analyses,  visual 
comparison  is  difficult  and  it  is  better  if  a  comparison  is  made  by  subtracting  calculations  at 
given  Doppler  velocity  values  along  x-axis.  The  result  of  doing  this  is  presented  in  Figure  4.7, 
which  indicates  an  absolute  difference  of  less  than  10'"^  for  all  Doppler  shifts. 
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Figure  4.7.  Difference  Plot  of  Zero  Time-Delay  Ambiguity  Surfaces  Generated  Using 
Independent  and  Dependent  Subcarrier  Induced  Doppler 


By  evaluating  both  plots,  one  is  inclined  to  use  the  more  efficient  IFFT  approach.  Even  if  this 
approach  does  not  give  a  perfectly  correct  output,  it  provides  a  much  faster  implementation  for 
generating  the  DAB  ambiguity  function.  This  speed  of  the  IFFT-based  algorithm  permits  faster, 
reliable  analysis  of  the  DAB  ambiguity  function  without  losing  considerable  information. 
However,  for  the  analysis  in  this  research  involving  filtering  the  reference  signal  from  the  target 
channel  of  a  PCL  receiver,  as  presented  in  a  later  section,  it  is  essential  to  model  the  signal 
completely  correct.  The  filter  designs  considered  are  a  function  of  the  signal  model  used  and 
require  a  perfect  (realistic)  signal  model. 
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4.3  DAB  waveform  performance  evaluation 


In  Chapter  3,  the  ambiguity  function  was  introduced  as  a  means  to  analyze  theoretical 
performance  of  a  waveform  of  opportunity.  The  width  of  the  ambiguity  function  peak  response 
is  a  measure  of  achievable  resolution  capability  and  accuracy  for  the  waveform  in  both  time 
delay  and  Doppler  measurements.  The  height  of  the  peak  response  is  a  measure  of  waveform 
detection  capability.  As  mentioned,  one  would  like  to  have  an  impulse-like  peak  response  having 
a  reasonable  height  above  the  noise  floor.  This  would  provide  reasonable  target  detection 
capability  with  good  resolution  in  both  time-delay  and  Doppler  frequency  measurements. 
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Figure  4.8.  Ambiguity  Surface  of  the  DAB  Waveform  (15  Subcarriers  and  10  Symbol  Coherent 

Integration  Time) 
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In  Figure  4.8,  the  DAB  waveform  performanee  is  characterized  through  the  calculation  of 
the  ambiguity  function.  In  this  case,  the  DAB  waveform  consisted  of  15  subcarriers  per  symbol 
and  a  10  symbol  coherent  integration  time.  In  this  normalized  plot,  the  surface  structure  near  the 
origin  (zero  Doppler  and  time-delay)  determines  waveform  capabilities.  The  narrower  the  peak 
response  and  neighboring  structure,  the  better  the  resolution  is  in  both  range  and  velocity.  More 
detailed  information  is  provided  in  Figure  4.9  and  Figure  4.10. 
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Figure  4.9.  Zero  Doppler  Ambiguity  Surface  Cut  of  the  DAB  Waveform  (15  Subcarriers  and  10 

Symbol  Coherent  Integration  Time) 

Figure  4.9  shows  the  zero  frequency  cut  of  the  DAB  waveform  ambiguity  surface.  In  this 
normalized  dB  plot,  the  maximum  appears  at  the  origin  as  expected.  Range  sidelobe  levels 
average  around  10.0  dB  below  the  maximum.  By  examining  the  plot,  one  can  get  an  idea  of  the 
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waveform  range  resolution  performance  as  well.  The  mainlobe  at  the  origin  covers  a  time 
interval  of  about  0.025  ms.  This  corresponds  to  about  3.75  km  of  monostatic  range-resolution 
which  may  be  considered  as  unsatisfactory  performance.  However,  it  should  be  remembered 
that  waveform  considered  for  this  illustration  is  highly  reduced  (lower  number  of  carriers  and 
symbols)  from  the  actual  DAB  waveform.  In  this  case,  one  DAB  symbol  consists  of  15 
subcarriers  whereas  the  actual  DAB  Mode  I  waveform  has  1536  subcarriers  per  symbol  as 
explained  in  Chapter  2.  It  should  also  be  noted  here  that  the  number  of  DAB  subcarriers  per 
symbol,  as  well  as,  the  number  of  symbols  per  coherent  integration  time  both  affect  the  radar 
performance  of  the  PCL  waveform.  The  relationship  between  mainlobe  (ML)  and  sidelobe 
levels  (SLL)  of  a  DAB  waveform,  having  M  subcarriers  and  N  symbols  per  coherent  integration 
time,  is  given  by  Equation  4.2  [53]. 

SLLdab  =ML-10*logl0(M*N) 
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Figure  4.10.  Zero  Time-Delay  Ambiguity  Surfaee  Cut  of  the  DAB  Waveform  (15  Subearriers 

and  10  Symbol  Coherent  Integration  Time) 

In  Figure  4.10,  the  zero  time-delay  eut  of  the  DAB  ambiguity  surfaee  of  Figure  4.8  is 
shown.  In  normalized  dB,  the  figure  shows  the  ambiguity  funetion  along  the  Doppler  frequeney 
axis  with  zero  time-delay.  As  illustrated,  the  main  Doppler  frequeney  sidelobes  are  about 
10.0  dB  below  the  peak  response  at  origin.  The  plot  also  provides  an  indieation  of  the  Doppler 
resolution  eapability  of  the  waveform.  A  Doppler  resolution  of  about  50  m/s  is  indieated  by  the 
data.  As  noted  previously  in  range  resolution  diseussions,  the  example  DAB  waveform  used 
here  for  illustration  is  sub-optimal  and  better  resolution  in  both  time  delay  and  Doppler  shift 
measurements  are  expeeted  from  a  waveform  with  an  inereased  number  of  subearriers  per 
symbol  and  symbols  per  eoherent  integration  time. 
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To  better  illustrate  DAB  performance,  more  simulated  data  was  generated  using  an 
increased  number  of  subcarriers  per  symbol  (75)  and  more  symbols  (50)  per  coherent  integration 
time.  Figure  4.11  is  the  zero-Doppler  frequency  cut  of  the  resultant  ambiguity  surface.  As 
mentioned  earlier,  this  ambiguity  surface  cut  represents  the  autocorrelation  of  selected 
waveform. 

dB 
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Figure  4.11.  Zero  Doppler  Ambiguity  Surface  Cut  of  the  DAB  Waveform  (75  Subcarriers  and 

50  Symbol  Coherent  Integration  Time) 

As  expected,  the  performance  here  is  better  than  the  previous  example  where  both  the 
number  of  subcarriers  and  symbols  per  coherent  integration  time  were  lower.  From  the 
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Figure  4.11,  the  sidelobe  levels  are  approximately  15  to  20  dB  below  the  peak  response  at  the 
origin  -  approximately  a  5.0  to  10.0  dB  improvement  over  the  previous  case. 
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Figure  4.12.  Zero  Time -Delay  Cut  of  DAB  Waveform  Ambiguity  Surface 
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Figure  4.12  is  the  zero  time-delay  cut  of  the  DAB  ambiguity  surface  for  the  new  DAB 
waveform.  This  axis  could  be  replaced  with  Doppler  shift  frequency  axis  for  the  case  where  the 
Doppler  shifted  signal  is  calculated  as  having  a  constant  Doppler  frequency  on  every  subcarrier. 
The  velocity  resolution  is  now  approximately  25  m/s  and  the  sidelobe  levels  are  around  15  to 
20  dB  lower  than  the  peak  at  the  maximum  -  approximately  a  5.0  to  10.0  dB  improvement  over 
the  previous  case. 
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The  zero-cut  ambiguity  surface  performance  results  presented  here  for  the  DAB 
waveform  closely  matches  calculations  available  from  the  literature  [18].  In  fact,  there  are  only 
a  few  open  literature  sources  for  comparison  of  results. 

The  DAB  waveform  performance,  as  supported  by  previous  data  plots  and  calculations 
above,  is  encouraging.  It  is  anticipated  that  the  accuracy  and  resolution  in  both  time  and 
frequency  measurements  will  be  quite  good  for  the  real  DAB  waveform,  given  it  consists  of 
1536  subcarriers  per  symbol.  In  practice,  there  will  be  certain  challenges  to  deal  with  in  order  to 
fully  exploit  benefits  of  the  DAB  waveform. 

4.4  DAB  detection  surfaces 

A  receiver  dedicated  to  exploiting  PCL  DAB  signals  may  use  a  matched  filter  for  signal 
processing  and  target  detection.  In  this  particular  signal-processing  scheme,  at  least  two  receiver 
channels  are  necessary.  One  channel  is  dedicated  to  receiving  the  reference  waveform,  defined 
as  the  direct  break-through  of  a  transmitted  signal  from  the  closest  line-of-sight  (LOS) 
transmitter  to  the  PCL  receiver.  The  second  receiver  channel  is  target  channel,  i.e.,  the  antenna 
of  this  receiver  channel  is  specifically  directed  towards  a  surveillance  region  to  receive  target 
returns  (echoes)  of  the  reference  waveform  returning  from  targets.  The  target  returns  are  simply 
time-delayed,  Doppler-shifted  versions  of  the  reference  waveform.  If  the  reference  waveform 
can  be  written  as  a  closed-form  mathematical  expression,  then  a  target  response  can  be  easily 
created  by  simply  incorporating  desired  time  delay  and  frequency  shift  to  resemble  a  target 
return. 

The  above-described  signal-processing  scheme  is  similar  to  the  process  implemented  for 
calculating  the  ambiguity  function,  as  explained  earlier.  The  ambiguity  surfaces  calculated 
earlier  represented  the  correlation  of  a  reference  DAB  waveform  with  a  time  delayed  and 
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Doppler  shifted  versions  of  itself.  Therefore,  the  peak  response  appeared  at  the  origin 
eorresponding  to  no  time  delay  and  no  Doppler  shift.  In  this  new  implementation,  the  reference 
waveform  is  both  time  delayed  and  frequency  shifted  before  being  correlated  with  the  target 
channel  output. 

In  general,  the  target  channel  output  consists  of  the  reference  signal  having  unknown  time 
delay  and  Doppler  frequency  shift,  as  induced  by  the  target  range  and  velocity  parameters.  In 
this  case,  the  peak  in  the  ambiguity  surface  appears  at  a  location  (time  delay  and  Doppler) 
corresponding  to  the  reference  waveform  having  the  same  phase  and  frequency  characteristics  as 
the  target  signal.  Therefore,  by  tracking  the  time-frequency  shifts  induced  on  reference 
waveform,  one  is  able  to  determine  the  corresponding  target  parameters  affecting  the  signal 
received  in  the  receiver  target  channel. 
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Figure  4.13.  Target  Channel  Detection  Surface  -  Target  Present  (30  Subcarriers  and  20  Symbol 

coherent  Integration  Time) 

Figure  4.13  is  a  target  channel  detection  surface  generated  to  illustrate  the  approach 
explained  above.  In  this  case,  the  DAB  waveform  consists  of  30  subcarriers  per  symbol  and  20 
symbols  per  coherent  integration  time.  The  simulated  target  channel  signal  in  this  example  only 
consists  of  a  target  response.  The  target  is  assumed  to  have  a  relative  velocity  of  200  m/s  in  a 
monostatic  sense,  i.e.,  the  target’s  radial  velocity  is  assumed  to  be  200  m/s  with  respect  to 
transmitter/receiver  site.  As  indicated,  the  surface  in  Figure  4.13  is  technically  not  an  ambiguity 
surface  but  rather  a  detection  surface,  i.e.,  a  surface  generated  using  techniques  similar  to 
ambiguity  function  calculations.  The  detection  surface  provides  theoretical  insight  into  how  well 
a  particular  waveform  may  perform  in  the  detection  and  estimation  of  targets.  In  this  case,  the 
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surface  is  generated  by  shifting  a  reference  signal  (in  both  time  and  frequeney)  and  eorrelating 
the  resultant  waveform  with  the  received  target  channel  signal.  The  relative  maximum(s)  of  the 
eorrelation  response  provides  a  measure  of  detection  and  estimation  performance  and  their 
relative  locations  are  closely  tied  to  the  target  range  and  veloeity  parameters. 
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Figure  4.14.  Target  Channel  Detection  Surface  -  Target  Present  (30  Subcarriers  and  20 

Symbol  Coherent  Integration  Time) 


For  elarity,  Figure  4.14  is  provided  and  represents  the  same  information  as  Figure  4.13 
but  from  a  different  visual  aspect  angle.  As  shown,  the  detection  surface  peak  occurs  at  the 
target’s  relative  velocity  location.  Therefore,  thresholding  eould  be  applied  to  deelare  targets 
present  -  hence,  the  terminology  detection  surface  is  applied. 
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In  practice,  the  received  target  channel  signal  does  not  only  include  target  returns  but  also 
noise.  The  model  chosen  to  reflect  the  ehannel  noise  contribution  is  Gaussian,  although,  the 
particular  choice  of  this  model  does  not  influence  research  results. 

Another  signal  component  present  in  the  received  target  channel  is  the  reference  signal 
itself  Therefore,  even  if  it  is  not  intended  to  have  the  reference  signal  at  the  target  channel 
output,  practieally  it  is  not  possible  to  isolate  the  target  channel  from  the  reference.  The  effeets 
of  this  component,  which  could  be  defined  as  the  leakage  of  the  reference  signal  into  the  target 
ehannel,  highly  degrade  the  performance  of  the  system.  Theoretieal  analysis  of  this  degradation 
is  implemented  in  Chapter  3. 

The  DAB  transmitter  architecture  is  built  in  aecordance  with  the  needs  of  Single  Frequency 
Networks  (SFN).  As  explained  earlier,  in  the  commercial  DAB  architecture  the  transmitters  are 
a  certain  distance  from  each  other  as  dietated  by  the  required  guard  interval  of  DAB  waveform 
(see  Chapter  2).  For  a  PCL  receiver,  this  means  that  in  the  received  signal,  along  with  the 
reference  signal  from  the  main  transmitter  closest  to  the  reeeiver,  there  are  time-delayed  replieas 
of  the  reference  signal  due  to  other  transmitters  in  the  vicinity.  These  replica  signals  are  lower  in 
power  since  they  are  from  transmitters  located  further  away  from  the  receiver.  For  this  research, 
these  replica  components  are  treated  as  nothing  more  than  a  multipath  reflection  of  the  reference 
signal  as  received  from  eaeh  distant  transmitter  in  the  SFN  structure. 

As  a  result  of  the  diseussion  above,  four  components  of  the  PCL  receiver  target  channel  can 
be  distinguished: 

1 .  The  return  of  the  transmitted  DAB  signal  from  the  target,  (this  eomponent  can  be 
thought  as  Doppler  and  time  (phase)  shifted  version  of  the  reference  signal). 
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2.  The  reference  signal  itself,  because  in  practice,  the  target  channel  cannot  be 
totally  isolated  from  the  direct  illumination  of  the  DAB  transmitter, 

3.  Multipath  response  of  the  reference  signal  due  to  multiple  transmitters 
simultaneously  transmitting  the  same  data  in  SFN  of  DAB.  The  delay  of  the 
multipath  signal  is  a  function  of  the  differential  distance  between  the  main 
transmitter  to  the  receiver  and  the  multipath  source  to  the  receiver, 

4.  The  channel  noise  component. 

Considering  the  received  signal  in  accordance  with  the  above  model,  significantly  affects 
the  detection  performance  of  a  DAB  signal  in  a  PCL  system.  Figure  4.15  represents  a  “DAB 
Detection  Surface  Plof  ’  for  a  received  signal  containing  all  four  components  discussed  above. 


Figure  4.15.  Representative  DAB  Detection  Surface  Plot 
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In  Figure  4.15,  the  reeeived  signal  eonsists  of  four  eomponents:  the  referenee  signal, 
multipath  component,  target  return  and  the  noise  component.  The  reference  waveform  consists  of 
30  subcarriers  and  20  symbols.  The  received  signal  carrier-to-Noise  Ratio  (CNR),  i.e.,  a  power 
ratio  of  reference  signal  and  noise,  is  30.0  dB.  The  multipath  component  power  level  at  the 
receiver  is  15.0  dB  lower  than  reference  signal.  The  target  return  power  level  in  the  received 
signal  is  30.0  dB  lower  than  the  carrier.  This  gives  an  SNR  value  of  around  0.0  dB  assuming  the 
target  return  component  is  the  “signal  of  interest”.  The  target’s  monostatic  equivalent  relative 
velocity  is  140  m/s. 

Figure  4.15  reveals  a  peak  response  on  the  surface  at  the  origin.  This  peak  response  is 
the  result  of  the  reference  and  multipath  component.  The  effect  of  this  high  correlation  is  not 
only  at  the  origin.  As  seen  in  the  ambiguity  function  analysis,  high  correlation  of  the  reference 
signal  with  itself  at  the  origin  extends  to  other  time-delay  and  frequency-shift  bins.  As  a  result, 
the  target  response  at  140  m/s  is  completely  “hidden”  by  the  reference  component  response  of 
the  received  signal  and  clearly  this  phenomenon  hinders  target  detection  for  a  CFAR  receiver. 
Therefore,  for  any  successful  matched  filter  detection  scheme,  where  time-delay  and  Doppler 
shift  of  a  target  signal  is  to  be  measured,  it  is  vital  to  eliminate  the  direct  path  signal  in  the  target 
channel. 

4.5  Direct  signal  filtering 

Several  approaches  may  be  taken  to  address  the  direct  signal  filtering  problem.  Different 
approaches  may  be  taken  in  different  domains.  The  first  two  domains  to  look  in  for  possible 
solutions  are  the  frequency  domain  and  the  time  domain.  Frequency  domain  solutions  attempt  to 
exploit  the  DAB  waveform  structure  in  frequency  domain  whereas  time  domain  solutions 
attempt  to  exploit  time  domain  properties  of  the  DAB  waveform. 
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An  initial  look  into  frequency  domain  techniques  did  not  seem  promising.  First,  the  DAB 
signal  has  a  relatively  large  bandwidth  of  1.5  MHz.  Any  filter  designed  in  frequency  domain 
would  need  to  differentiate  between  the  reference  signal  and  the  target  component  of  the 
received  signal  in  the  frequency  domain.  The  target  return  in  frequency  domain  basically 
appears  as  a  shifted  version  of  reference  signal,  with  the  amount  of  Doppler  frequency  shift  is 
associated  with  target  relative  velocity. 


Figure  4.16.  Frequency  Domain  Response  of  DAB  Direct  Signal  and  Doppler  Shifted 

Target  Return 


One  can  infer  from  Figure  4.16  that  for  reasonable  Doppler  frequencies,  i.e.,  those 
generated  from  velocities  of  typical  airborne  platforms,  it  is  not  possible  to  use  narrowband- 
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filtering  techniques  of  traditional  radar  signal  processing  to  eliminate  the  reference  signal 
response.  Traditional  solutions  including  frequency  domain  HR  and  FIR  filtering  seem 
impractical  for  eliminating  the  reference  signal  while  passing  the  target  return  unaltered. 

If  working  in  the  time  domain,  more  interesting  approaches  can  be  considered.  One  idea 
involves  using  intelligent  signal  processing  schemes  from  the  world  of  Artificial  Neural 
Networks  (ANN)  [24]. 

Another  unique  solution  was  implemented  as  part  of  this  research  and  comes  from  the 
basic  idea  of  a  matrix  projection  operation.  Detailed  background  material  about  this  projection 
operation  is  presented  in  Chapter  3.  Basically,  the  idea  is  to  “project”  the  received  signal  into  an 
alternate  space  that  is  orthogonal  to  the  direct  path  signal.  Ideally,  this  projection  eliminates  the 
direct  path  component  from  the  projected  signal  while  having  minimal  affect  on  the  desired 
target  return.  Target  detection  is  then  implemented  using  the  resultant  signal  from  the  projection 
operation  of  the  received  signal.  There  are  several  questions  that  have  to  be  answered  in  this 
signal-processing  algorithm: 

1.  Does  the  projection  operator  completely  eliminate  the  direct  path  signal  consisting  of 
a  reference  component  together  with  multipath  component? 

2.  Does  the  projection  preserve  the  desired  target  return  component? 

3.  How  can  the  projection  operation  be  implemented  for  the  DAB  waveform? 

What  one  wants  is  a  design  that  completely  eliminates  the  reference  and  multipath 
component,  where  the  multipath  component  is  simply  a  time-delayed  version  of  the  reference 
component.  It  is  also  desirable  for  the  target  return  component  of  the  received  waveform  to  pass 
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through  the  filtering  unchanged.  Of  course,  these  are  ideal  considerations  and  practically  any 
realizable  filter  implementation  will  not  be  able  to  completely  eliminate  the  direct-path  signal. 
As  with  all  filter  designs,  the  projection  operation  technique  will  probably  attenuate  target 
returns  as  well. 

Considering  the  discussion  on  projection  operation  effects,  a  new  set  of  goals  is 
introduced.  Namely,  one  wants  the  direct-path  signal  to  be  sufficiently  attenuated  such  that  its 
response  on  the  detection  surface  does  not  hinder  target  detection,  i.e.,  it  does  not  have  a 
response  that  “hides”  the  desired  response  at  time  delay  and  Doppler  shift  values  consistent  with 
the  target’s  parameters.  This  effect  was  illustrated  in  Figure  4.15  and  subsequently  discussed. 

Another  important  criterion  for  evaluating  practical  filter  performance  is  the  evaluation  of 
Signal-to-Noise  Ratios  (SNR)  before  and  after  filtering.  For  this  work,  the  “signal  of  interest”  is 
the  target  return  component  of  the  received  signal.  Although  there  will  be  signal  attenuation 
during  filtering,  the  resultant  SNR  values  after  filtering  should  be  at  reasonable  levels  to  enable 
target  detection  and  parametric  estimation.  This  brings  another  consideration  into  play,  namely, 
what  is  the  projection  operation  effect  on  the  noise  component  of  the  received  signal? 

Details  of  the  proposed  time  domain  projection  filter  are  included  in  Chapter  3,  where 
information  about  the  exploitation  of  DAB  signal  structure  to  form  the  projection  operation  is 
explained.  Simulation  results  for  the  proposed  projection  filter  approach  are  presented  next. 

4.6  Simulations 

The  simulation  results  presented  in  this  section  are  for  two  different  scenarios.  The  first 
scenario  is  an  ideal  scenario  where  the  phase  reference  symbol  transmitted  in  every  frame  is 
demodulated  under  ideal  conditions,  i.e.,  the  phase  reference  for  every  subcarrier,  as  used  by  the 
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receiver  to  demodulate  and  estimate  bits,  is  free  of  error.  The  second  scenario  is  more  realistic 
and  based  on  imperfect  extraction  and  recovery  of  the  phase  reference  information. 

4.6.1  Ideal  case:  Perfect  phase  references 

Having  perfect  phase  references  in  the  ideal  scenario  provides  a  perfect  base  for  the  filter 
implementation.  A  perfect  base,  as  discussed  in  Chapter  3,  is  a  signal  consisting  of  the 
superposition  of  the  reference  signal  and  the  multipath  component.  The  reference  and  multipath 
components  both  consist  of  1536  subcarriers  each.  As  shown,  this  subcarriers  may  be 
represented  via  complex  phasor  representation  and  the  complex  phasors  are  added  vectorially  to 
create  the  sum.  Each  subcarrier  has  a  certain  phase  code  in  conjunction  with  the  data  that  is 
being  modulated.  For  phasor  addition  of  subcarriers,  the  phase  code  and  amplitude  of  the 
reference  subcarrier  represents  one  phasor  and  the  phase  (most  likely  different  from  reference 
component  subcarrier)  and  an  amplitude  (smaller  due  to  longer  propagation  path  length)  of  the 
multipath  component  represent  the  other.  A  perfect  base  consists  of  the  superposition  of  both 
phasor  components  for  each  subcarrier.  By  developing  a  filter  design  using  a  perfect  base, 
meaning  the  received  signal  is  projected  into  a  space  orthogonal  to  the  perfect  base  according  to 
Chapter  3,  the  filter  is  able  to  completely  remove  the  direct-path  component  from  the  received 
signal. 

Under  conditions  of  the  ideal  scenario,  we  expect  to  completely  remove  the  effects  of  the 
direct-path  component  from  the  detection  surface.  After  filtering,  the  resultant  signal  is 
correlated  with  a  time  delayed  and  frequency  shifted  version  of  the  reference  signal  to  form  the 
detection  surface.  Ideally,  the  remaining  response  on  the  detection  surface  is  due  solely  to  the 
target  return  and  its  parameters  can  be  estimated. 
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Figure  4.17.  Post-Filtered  Detection  Surface  Showing  Detectable  Target  Response  Ideal 

Scenario  -  Perfect  Symbol  Phase  Reference 


Figure  4.17  is  the  post-fdtered  detection  surface  for  ideal  conditions,  i.e.,  symbol  phase 
references  are  perfectly  known.  All  radar  and  signal  parameters  for  this  illustration  are 
consistent  with  those  used  for  creating  Figure  4.15.  Clearly,  the  direct  path  response  at  the  origin 
is  almost  totally  removed.  This  resultant  suppression  at  the  origin  is  due  to  filtering  the  received 
signal,  where  the  ideal  filter  effectively  removes  both  the  reference  and  multipath  component  of 
the  received  signal.  For  clarity.  Figure  4.18  is  provided  and  represents  the  same  information  as 
Figure  4.17  but  from  a  different  visual  aspect  angle.  As  shown,  the  detection  surface  peak 
occurs  at  the  target’s  relative  velocity  location.  Therefore,  thresholding  could  be  applied  to 
declare  targets  present  -  hence,  the  terminology  detection  surface  is  applied. 
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Figure  4.18.  Post-Filtered  Detection  Surface  Showing  Detectable  Target  Response  Ideal 
Scenario  -  Perfect  Symbol  Phase  Reference 

4.6.2  Realistic  conditions;  Imperfect  phase  references 

After  seeing  the  effectiveness  of  a  filter  designed  under  ideal  conditions,  more  realistic 
conditions  are  simulated  for  the  second  scenario,  i.e.,  the  phase  reference  symbol  cannot  be 
recovered  error-free.  In  practical  applications,  this  is  always  true  because  the  received  phase 
symbol  consists  of  four  components:  reference  signal,  multipath  signal,  target  return  and  noise. 
Therefore,  the  thresholds  calculated  from  the  phase  reference  symbol  will  not  be  perfect.  The 
perfect  base  signal,  which  was  perfectly  generated  from  the  superposition  of  reference  and 
multipath  components,  is  no  longer  available  under  scenario  two.  This  is  true  because  each 
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received  subcarrier  is  now  the  superposition  of  a  reference  signal,  multipath  components,  a  target 


return  (assuming  a  target  is  present)  and  a  channel  noise  component.  This  is  shown  in  Chapter  2. 


The  filter  design  implemented  using  this  new  imperfect  base  signal  will  not  be  able  to 


totally  remove  the  direct-path  component.  The  transformation  space  defined  by  the  orthogonal 


projection  of  the  received  signal  is  now  orthogonal  to  the  imperfect  base  signal  comprised  of 


components  other  than  the  direct-path.  This  effectively  causes  a  leakage  of  the  reference  and 


multipath  components  into  the  projected  signal,  as  well  as,  more  significant  attenuation  of  the 


desired  target  return. 
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Figure  4.19.  Pre-Filtered  Detection  Surface  Showing  Undetectable  Target  Response 
Non-Ideal  Scenario  -  Imperfect  Symbol  Phase  Reference 
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Figure  4.19  is  the  pre-filtered  detection  surface  for  non-ideal  conditions,  i.e.,  symbol 
phase  references  are  not  perfectly  known.  As  before,  the  received  signal  has  a  direct  signal 
reference  component  that  is  30.0  dB  above  both  the  target  return  and  the  noise  component,  and  it 
is  10.0  dB  higher  than  the  multipath  component.  Because  of  direct  path  dominance  in  the 
received  signal,  the  only  discernable  peak  appears  near  the  origin,  just  like  an  ambiguity  surface 
plot.  For  clarity.  Figure  4.20  is  provided  and  represents  the  same  information  as  Figure  4.19  but 
from  a  different  visual  aspect  angle.  As  shown,  the  unfiltered  response  near  the  origin 
completely  dominates  and  the  target  response  is  undetectable. 
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Velocity 
in  m/s 


Figure  4.20.  Pre-Filtered  Detection  Surface  Showing  Undetectable  Target  Response 
Non-Ideal  Scenario  -  Imperfect  Symbol  Phase  Reference  (side  view) 
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Figure  4.21,  is  the  plot  of  the  post-filtered  deteetion  surface.  A  better  plot  showing  the 
effect  of  filtering  is  Figure  4.22.  Care  must  be  taken  in  analyzing  this  example,  the  direct  path 
response  on  the  detection  surface  is  not  completely  removed. 
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Figure  4.21.  Post-Filtered  Detection  Surface  Showing  Detectable  Target  Response  Non- 
Ideal  Scenario  -  Imperfect  Symbol  Phase  Reference 


One  encouraging  fact  that  shows  the  efficiency  of  this  filtering  approach  is  that 
the  target  response  is  now  visible.  The  target  detection  performance  increase  is  a  result 
of  filtering  the  received  signal  direct-path  component.  Thus,  a  qualitative  performance 
analysis  is  necessary  to  see  the  filter  design  effectiveness. 

For  this,  an  example  situation  is  considered  where  the  received  signal  consists  of 
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the  reference  signal,  a  multipath  component  and  noise.  For  this  analysis,  the  DAB 
waveform  of  opportunity  consists  of  30  subcarriers  per  symbol  and  20  symbols  per 
coherent  integration  time. 
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Figure  4.22.  Post-Filtered  Detection  Surface  Showing  Detectable  Target  Response  Non- 
Ideal  Scenario  -  Imperfect  Symbol  Phase  Reference  (side  view) 


Figure  4.23  is  the  spectral  response  of  a  DAB  signal.  The  solid  line  represents  the 
spectral  estimate  obtained  from  an  averaged  periodogram  [25].  The  plot  is  normalized  and 
shows  amplitudes  as  a  function  of  frequency  in  dB.  The  dashed  line  represents  the  spectral 
estimate  after  filtering.  Clearly,  the  amount  of  reference  signal  power  passing  through  the  filter 
is  sharply  reduced  (approximately  22  dB  reduction  in  peak  response).  This  reduction  in 
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reference  signal  power  is  the  phenomenon  that  enables  the  target  return  component  of  the 
received  signal  to  be  detected  using  the  detection  surface  technique  in  this  research. 
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Figure  4.23.  Spectral  Filter  Responses  for  DAB  Signal 
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4.7  Performance  evaluation  of  direct-path  filter 

The  results  presented  in  Figure  4.23  should  be  interpreted  with  caution.  The  proposed 
filter  for  this  research  projects  the  received  signal  into  a  space  orthogonal  ideally  to  the  direct- 
path  component  of  the  received  signal.  The  underlying  assumption  is  that  this  projection  totally 
eliminates  the  direct-path  response.  However,  as  explained  above,  the  algorithm  does  not  create 
a  perfect  basis  for  ideal  projection,  due  to  the  target  and  noise  components  of  the  received  signal. 
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Also,  the  target  eomponent  itself  is  not  perfeetly  orthogonal  to  the  direet-path.  Thus,  the  target 
eomponent  will  be  attenuated  by  this  projeetion  operation,  as  well  as,  the  noise  eomponent. 
Therefore,  a  qualitative  analysis  of  various  attenuation  levels  for  signal  eomponents  is  neeessary 
and  follows. 

In  Figure  4.20,  a  pre-filtered  deteetion  surface  is  plotted  for  the  received  signal.  The  peak 
response  level  occurs  at  the  origin  and  is  approximately  31.47  dB.  When  considering  the 
detection  surface  response  at  the  origin  after  filtering,  as  in  Figure  4.22,  the  level  is  about 
12.55  dB.  Therefore,  a  reduction  of  18.92  dB  can  be  associated  with  the  direct  path  response  of 
the  detection  surface. 

Similar  analysis  can  be  done  for  a  target  response  on  the  detection  surface.  Before 
filtering,  the  target  response  was  22.12  dB.  After  filtering,  the  detection  surface  is  similar  to 
Figure  4.22  and  the  target  response  is  approximately  14.31  dB.  The  target  response  is  reduced 
by  approximately  7.81  dB. 

Similarly,  the  noise  floor  is  attenuated  by  the  filtering  process  and  causes  a  change  in 
SNR  after  filtering.  As  a  measure  of  noise  floor  attenuation,  every  pre-filtered  detection  surface 
component  is  added  and  compared  to  the  sum  of  detection  response  components  after  filtering. 
When  done,  the  attenuation  of  the  noise  floor  attenuation  is  around  1 1.0  dB.  This  simple  analysis 
technique  provides  an  idea  of  potential  noise  floor  reduction,  although  a  rigorous  proof  of  the 
same  is  not  provided. 

The  numbers  calculated  highlight  the  effectiveness  of  the  filter.  Although  all  received 
signal  components  are  attenuated,  including  the  desired  target  return,  the  interfering  direct  path 
component  is  attenuated  much  more  significantly.  Effectively,  the  higher  direct  path  attenuation 
results  into  an  11.11  dB  accentuation  of  the  target  response,  relative  to  the  direct  path  response. 
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on  the  detection  surface.  Also,  the  detection  surface  noise  floor  is  attenuated  more  than  target 
response. 

4.8  Summary 

In  this  chapter  the  proposed  direct  path  filtering  approach  introduced  in  Chapter  3  was 
simulated  and  its  performance  characterized.  It  was  found  to  effectively  eliminate  (or 
significantly  reduce)  the  direct-path  signal  in  the  target  channel  of  a  PCL  receiver  containing 
direct-path,  target  return  and  noise  components. 

The  chapter  started  with  the  performance  evaluation  of  DAB  waveform  as  a  PCL 
waveform  of  opportunity.  It  has  been  found  that  with  DAB  has  encouraging  ambiguity  function 
properties.  Time-delay  and  Doppler  shift  measurement  accuracy  and  resolution  can  be  dealt 
independently  and  both  can  be  improved  simultaneously.  Range  and  frequency  sidelobe  levels 
are  found  to  be  a  function  of  the  number  of  subcarriers  and  the  number  of  symbols  per  coherent 
integration  time.  Equation  4.3  has  been  shown  to  be  consistent  with  theoretical  calculations 
implemented  for  this  research. 

SLLdab  =ML-10*logl0(M*N) 

In  Equation  4.3,  SEE  is  the  DAB  sidelobe  level  in  dB,  ML  is  the  mainlobe  level  in  dB,  M 
is  the  number  of  subcarriers  per  symbol  and  N  the  number  of  symbols  per  integration  time. 

A  serious  problem  in  a  practical  DAB  PCL  system  implementation  is  identified,  i.e.,  the 
existence  of  a  strong  direct-path  component  in  the  PCL  receiver  target  channel.  A  unique  time 
domain  filtering  technique  has  been  proposed  to  eliminate  this  (Chapter  3)  and  its  performance 
evaluated  in  Chapter  4. 
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The  filter  design  performance  is  evaluated  for  two  different  scenarios,  an  ideal  scenario 
based  on  perfect  phase  reference  estimation  and  a  non-ideal  scenario  based  on  imperfect  phase 
reference  estimation.  It  was  determined  that,  although  the  proposed  filter  design  could  not 
totally  eliminate  the  direct  path  signal  from  the  target  channel,  it  did  significantly  attenuate  the 
direct-path  signal  such  that  the  relative  target  response  on  the  detection  surface  could  be 
detected. 

The  next  chapter  includes  some  important  conclusions  and  considerations  about  the  filter 
design  and  possible  future  work  to  expand  upon  this  research. 
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CHAPTER  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Summary 

This  research  focused  on  using  Digital  Audio  Broadcast  (DAB)  signals  within  a  Single 
Erequency  Network  (SEN)  architecture  for  Passive  Coherent  Location  (PCL)  applications  -  the 
DAB  signal  promises  excellent  results  as  a  radar  waveform  [56].  In  this  work,  the  advertised 
DAB  performance  as  a  waveform  of  opportunity  in  a  PCL  applications  is  verified. 

More  importantly,  a  serious  problem  of  practical  PCL  systems  using  matched  filters  to 
extract  the  time-delay  and  Doppler-shift  information  is  addressed.  The  problem  arises  from  the 
unwanted  direct-path  signal  component  in  the  receiver  target  channel.  Unless  completely 
eliminated  or  significantly  attenuated,  this  component  renders  target  detection  via  simple 
matched  filtering  impossible.  The  critical  effect  of  the  reference  signal  component  in  the  target 
channel  of  the  matched  filter  receiver  is  shown. 

To  deal  with  the  problem,  a  time-domain  filter  design  is  proposed  which  exploits  the 
DAB  orthogonal  signal  structure  to  eliminate  the  reference  component  from  the  target  channel 
signal.  The  filter  performance  is  evaluated  for  two  cases,  including: 

1 .  An  ideal  (unrealistic)  case,  where  the  phase  reference  symbol  for  every  frame  of 
DAB  signal  is  assumed  to  be  recovered  error-free.  This  case  provided  the  theoretical 
limit  of  filter  performance. 

2.  A  more  realistic  case  where  channel  frequency  response  affects  the  complete 
DAB  signal,  as  expected  in  a  realistic  system.  This  case  provided  the  practical 
capabilities  of  the  filter. 
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Results  for  the  ideal  case  showed  perfect  elimination  of  the  direct-path  component  from 
the  target  channel  signal,  resulting  in  robust  detection  of  target  returns  as  shown  through  creation 
of  detection  surfaces  before  and  after  filtering  the  received  signal  with  the  proposed  filter. 

For  case  two,  the  direct-path  component  was  not  totally  eliminated.  However,  results 
show  that  the  proposed  filter  design  did  attenuate  the  direct-path  component  sufficiently  such 
that  target  returns  were  detected  using  matched  filtering.  The  proposed  filter  design  is  thus 
declared  successful. 

5.2  Conclusions 

As  proven  in  this  research,  the  proposed  time-domain  filter  for  removing  the  direct-path 
signal  from  the  DAB-PCL  receiver  target  channel  is  effective  and  allows  target  detection  via 
matched  filtering.  Other  issues  encountered  throughout  the  research  are  addressed  here  in  the 
conclusion. 

The  proposed  projection  filtering  operation  effectively  removes  the  direct-path 
component  from  the  receiver  target  channel.  However,  the  projection  filter  operation  also 
attenuates  both  the  target  return  and  noise  components.  This  means  more  sensitive  receivers  may 
be  required  and  increased  sensitivity  may  present  a  challenge  for  current  systems.  Nevertheless, 
the  SNR  value,  a  ratio  of  target  return  signal  strength  to  noise  power,  did  not  decrease.  This 
means  target  detection  performance  will  not  be  negatively  impacted  provided  the  receiver 
sensitivity  does  not  become  a  major  issue  because  of  the  received  signal  component  attenuation. 

The  initial  impression  for  the  target  return  component  in  the  target  channel  is  that  one 
would  want  this  component  to  be  as  strong  as  possible,  relative  to  the  direct-path  component  in 
the  same  channel.  This  is  due  to  the  fact  that  increased  target-return  level  means  a  larger 
response  on  the  detection  surface  at  a  time-delay  and  Doppler  shift  value  corresponding  to  true 
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target  parameters.  However,  another  important  faetor  for  robust  target  deteetion  is  the  filter 
effieiency  in  eliminating  the  direet-path  component  from  the  target  channel  response.  Filter 
efficiency  is  a  function  of  direct-path  signal  strength. 

As  shown  in  Chapter  4,  the  filter  is  ideally  designed  for  the  case  when  the  phase 
reference  symbol  in  the  DAB  transmission  frame  is  perfectly  demodulated  resulting  in  error-free 
phase  thresholds  for  each  subcarrier  needed  for  the  filter  design.  The  filter  (detection) 
performance  of  the  system  degrades  when  noise  and  the  target  return  component,  is  added  to  the 
received  signal  as  expected  realistic  scenarios.  Therefore,  an  increased  target  return  component 
in  the  target  channel  does  not  necessarily  mean  better  system  detection  performance  because  it 
can  cause  less  optimal  filter  designs  resulting  in  less  efficient  filtering  of  the  direct-path  in  the 
target  channel. 

A  similar  argument  can  be  presented  with  regard  to  the  direct-path  signal  level,  where  the 
direct-path  signal  is  the  sum  of  the  reference  and  multipath  components.  Increasing  the  direct- 
path  signal  level  in  the  target  channel  does  not  necessarily  mean  decreasing  system  detection 
performance,  since  for  the  proposed  filtering  scheme  the  robustness  and  efficiency  increases 
when  dominant  direct-path  signal  levels  are  present. 

The  DAB  waveform  used  in  this  research  was  significantly  downscaled  from  the  actual 
waveform.  The  actual  DAB  Mode  I  signal  consists  of  1536  subcarriers.  The  highest  number  of 
subcarriers  used  in  this  research  is  75.  In  addition,  coherent  integration  time  was  also 
downscaled  to  include  a  manageable  number  of  symbols  per  integration  time.  Therefore,  the 
performance  of  the  real  DAB  waveform  is  expected  to  be  better  [18,19]. 
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5.3  Recommendations  for  Future  Research 


This  research  represents  one  step  in  enabling  the  use  of  a  DAB  waveform  of  opportunity 
in  a  passive  radar  application  using  a  matched  filter  receiver.  When  the  research  for  this  effort 
was  conducted  several  topics  emerged  for  possible  future  research  work  and  are  presented. 

The  first  and  most  related  future  study  could  include  using  the  proposed  projection  filter 
design  with  real-time  data.  The  collected  data  could  be  fed  into  a  simulated  PCL  receiver.  After 
the  received  signal  is  filtered,  realistic  detection  schemes  could  be  evaluated  to  determine 
performance  of  the  filtered  DAB  signal.  Closed-form  solutions  for  probability  of  false  alarm 

and  probability  of  detection  may  be  developed  and  a  threshold  for  reasonable  values 

calculated.  This  follow-on  work  could  be  very  helpful  in  showing  the  performance  of  the  filter 
design  with  real  data. 

The  proposed  filter  design  of  this  work  exploits  the  time  domain  structure  of  the  DAB 
signal.  The  DAB  time  domain  structure  may  be  exploited  in  other  ways.  One  effort  may 
concentrate  on  the  constant  structure  of  DAB  transmission  frame.  Namely,  every  DAB  frame 
starts  with  a  null  symbol  defining  the  start  of  the  frame.  The  null  symbol  is  followed  by  the 
phase  reference  signal  having  a  pre-determined  modulation  known  to  the  receiver.  Modulations 
for  phase  reference  symbols  remain  constant.  Thus  for  slowly  changing  channel  conditions,  a 
fixed  and  constant  symbol  is  present  within  every  frame  following  the  null  symbol. 

The  stable  features  of  the  null  and  phase  reference  symbols  may  be  seen  as  a  constant 
repeating  structure  occurring  within  every  frame.  This  suggests  there  may  be  potential  to  use  the 
constant  repeating  signal  structure  just  like  pulse  repetition  is  used  in  pulse  radar  applications. 
There  will  be  a  pulse  repetition  frequency  associated  with  the  waveform,  which  can  be  defined  as 
the  reciprocal  of  frame  time  since  the  constant  structure  repeats  once  in  every  frame.  With  frame 
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times  being  around  100  milliseeonds,  the  resultant  pulse  repetition  frequeney  would  be  very  low 
and  lead  to  very  ambiguous  Doppler  measurements.  However,  at  the  same  time  this  potentially 
means  very  unambiguous  range  measurement  eapability.  Coherent  integration  time  may  be 
ehosen  to  inelude  any  number  of  desired  pulses.  The  number  of  pulses  eoherently  proeessed  ean 
signifieantly  affeet  the  system’s  eapability  to  deteet  targets.  If  a  low  number  of  pulses  are  used, 
Pfa  and  eriteria  might  not  be  met.  However,  if  a  large  number  of  pulses  per  eoherent 

integration  time  is  used,  the  system  may  have  an  undesirable  slow  response  time.  These  ideas 
and  other  related  issues  eoneerning  pulsed  PCL  system  applieations  may  be  worthy  of  future 
study  efforts. 
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APPENDIX  A  -  Acronyms 


ANN 

Artificial  Neural  Networks 

ARM 

Anti-radiation  Missile 

BER 

Bit  Error  Rate 

GEAR 

Constant  False  Alarm  Rate 

COEDM 

Coded  Orthogonal  Frequeney  Division  Multiplex 

DAB 

Digital  Audio  Broadeast 

DBS 

Direet  Broadeast  Satellite 

DOA 

Direetion  of  Arrival 

DQPSK 

Differential  Quadrature  Phase  Shift  Keying 

ECM 

Eleetronie  Counter  Measures 

EMCON 

Emission  Control 

ESM 

Eleetronie  Support  Measures 

ETSI 

European  Teleeommunieations  Standards  Institute 

EDM 

Frequeney  Division  Multiplex 

EEC 

Forward  Error  Correetion 

FEE 

Fast  Fourier  Transform 

FIR 

Finite  Impulse  Response 

EM 

Frequeney  Modulation 

ICI 

Inter  Carrier  Interferenee 

IFF 

Identifieation  Friend  or  Foe 

IFFT 

Inverse  Fast  Fourier  Transform 

HR 

Infinite  Impulse  Response 
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ISI 

Inter  Symbol  Interference 

OFDM 

Orthogonal  Frequency  Division  Multiplex 

PCL 

Passive  Coherent  Location 

RF 

Radio  Frequency 

SFN 

Single  Frequency  Network 

SNR 

Signal  to  Noise  Ratio 

TDOA 

Time  Difference  of  Arrival 

VHF 

Very  High  Frequency 
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APPENDIX  B  -  Simulation  Code 


function[]=ftest (subi,c) ; 

g _ 

0 - 

-g  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 


%Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
%At  Air  Force  Institute  of  Technology 
%Dayton-OH  USA  2002 

^ ================================================= 

g _ 

0 - 

%master  simulation  code 


[vec  ref  recsig] =recvec (subi, c) ; 

f symbol=ref ( 1 : length (recsig) / c)  ; 
p_mat=reffaz (fsymbol) ; 
DEMOD=demode2 (p_mat) ; 

XF=subs (DEMOD) ; 
mysig=f iltre (XF, recsig)  ; 

figure; 

ambfun (recsig)  ; 
figure; 

ambfun (mysig) ; 
return 
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function [vec, ref, signal] =recvec (subi,  c)  ; 


g _ 

0 - 

g _ 

0 - 

%Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
%At  Air  Force  Institute  of  Technology 
%Dayton-OH  USA  2002 

g _ 

0 - 

g _ 

0 - 

%Received  signal  creator 
if  nargin~=2 

error (' please  use  testi (subi, c)  format'); 

end 

vec=sembol (subi, c) ; 

ch=input (' target  echo  wanted  ? (type  anything  if  no):'); 
load  sembolde; 
if  isempty (ch) ==1 

doo=input (' please  enter  the  doppler  velocity  required  for  dss:'); 
dsembol2 (subi, c, doo) ; 
load  sembolde; 

sdr=input (' please  enter  SDR (  signal  to  doppler  echo  ratio  ):'); 

ayar (dvec2, sdr) ; 

load  ayarde; 

dvec2=signal ; 

else 

dvec2=zeros (size (vec) )  ; 

end 

save  dopsigde  dvec2 

ch=input('do  you  want  multipath  ?  (  type  anyting  if  no) :'); 

if  isempty (ch) ==1 

mut=input ( ' give  a  multipath  between  0-0. 24 6ms  : ' ) ; %actually  it  becomes 
:between  0-50  percent 
coky=multip (vec,mut)  ; 

smr=input (' please  enter  SMR  (  signal  to  multipath  ratio  ) :'); 

ayar (coky, smr) ; 

load  ayarde; 

coky=signal ; 

else 

coky=zeros (size (vec) ) ; 

end 

ch=input('do  you  want  gaussian  noise  ?  (  type  anything  if  no  ):'); 

if  isempty (ch) ==1 

gurultu (vec) ; 

load  gurultude; 

else 

nn=zeros (size (vec) ) ; 

end 

signal=vec+dvec2+coky+nn ; 

ref=vec+coky; 

return 
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function [vec] =sembol (subi, c) ; 


g, 

0 

g, 

0 


%Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
%At  Air  Force  Institute  of  Technology 
%Dayton-OH  USA  2002 


g, 

0 

g, 

0 


%DAB 


symbol 


creator 


if  nargin~=2 

error ('use  format  sembol (subi, c)  ' )  ; 

end 


fc  =  200000000; 
totalt  =  .001; 
vec=  [  ]  ; 

f reqc= [fc:1000: (subi-l)*1000+fc]; 
f req=f reqc-f c+1000 ; 
f s=2*max (freq) ; 
sample_n=totalt*f s ; 

Fs=2^ (3+length (dec2bin (sample_n) ) ) ; 

Fc=median (freq) +fc; 

delta=totalt /Fs ; 

t= [ 0 : delta : totalt-delta]  ; 

rn=floor (rand (c, subi) *4) ; 

save  random  rn 

rnn=pi/2*rn; 

mo=exp ( j*rnn) ; 

coef= [ zeros (c, 1 )  mo  zeros (c, Fs-subi-1 )] ; 
for  p=l:c; 

icoef=real (ifft (coef  (p,  : ) ) ) ; 
nicoef=icoef /max (icoef ) ; 
vec=[vec  nicoef ] ; 
end 

save  parameter  c  subi  Fs  totalt  freq  t 

save  sembold 

save  sembolde  vec; 

return 
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function [ ] =dsembol2 (subi,  c,  dop)  ; 


g _ 

0 - 

g _ 

0 - 

%Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
%At  Air  Force  Institute  of  Technology 
%Dayton-OH  USA  2002 

g _ 

0 - 

g _ 

0 - 

%Doppler  shifted  signal  creator 
if  nargin~=3 

error('use  this  model : dsembol (subcarrier, symbol, doppler  in  m/s)'); 

end 

fc  =  200000000; 
totalt  =  .001; 
dvec2= [ ] ; 

f reqc= [ f c : 1000 :  ( subi-1 ) *1000  +  fc]  ; 
f req=f reqc-f c+1000  ; 
f s=2*max (freq) ; 
sample_n=totalt*f s ; 

Fs=2^ (3  +  length (dec2bin (sample_n) ) )  ; 

Fc=median (freq) +fc; 

load  random 
delta=totalt /Fs ; 

t= [ 0 : delta : tot alt-delta] ; %time  vector 
fd  =  2*dop*Fc/3/10^8; 

z  =  l; 

for  p=l:c; 

pt=t+ (p-1) *totalt; 

for  i=l : length ( freq) ; 

x(z,  :)=cos(2*pi*freq(i)*t+(pi/2)*rn(p,i) ) ; 

y(i,  :)=cos(2*pi*freq(i) *t+2*pi*fd*pt+ (pi/2) *rn (p, i) ) ; 

z=z+l ; 

end 

dopp_sigp=sum (y,  1 )  ; 

dopp_sigp=dopp_sigp/max (dopp_sigp)  ; 
dvec2=[dvec2  dopp_sigp] ; 

end 

save  sembolde  dvec2; 
save  nullspace  x; 
return 
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function [ ] =ayar (signal,  snr)  ; 


g _ 

0 - 

g _ 

0 - 

%Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
%At  Air  Force  Institute  of  Technology 
%Dayton-OH  USA  2002 

g _ 

0 - 

g _ 

0 - 

%Signal  to  noise  ratio  calculator 
if  nargin~=2 

error  ('use  the  format  ayar (signal, attenuation) ,  at.  in  dB ' )  ; 
end 

if  isempty (snr) ==1 
ratio=l ; 

else 

ratio=10^ (snr/10) ; 

end 

signal=signal/sqrt (ratio)  ; 
save  ayarde  signal; 
return 


no 


function [mulsig] =multip (sinyal,  delay)  ; 


%Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
%At  Air  Force  Institute  of  Technology 
%Dayton-OH  USA  2002 


g, 

0 

g, 

0 


%multipath 


signal 


creator 


if  nargin~=2 

error ('enter  two  parameters : multip (vector , multipath  in  ms)  :1ms 
corresponds  to  150  km  range  multipath  so  be  careful  \n  range  between  0 
to  .5ms  \n  (  if  you  have  an  arbitrary  signal  then  this  range  is  0-.5  of 
vectors  length'); 
end 

load  parameter; % 
kk=length (sinyal ) / c ; 
mulsig= [ ] ; 
for  i=0:c-l; 

kak=cokyol (sinyal ( l  +  i*kk :  ( 1+i ) *kk) , delay )  ; 
mulsig= [mulsig  kak]  ; 

end 

return 
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function [xf ] =subs (MM) ; 

g _ 

0 - 

g _ 

0 - 

%Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
%At  Air  Force  Institute  of  Technology 
%Dayton-OH  USA  2002 

g _ 

0 - 

g _ 

0 - 

%subcarrier  calculator 

load  parameter 

z  =  l; 

for  p=l:c; 

for  i=l : length ( freq) ; 

xf (z,  : )=cos (2*pi*freq(i) *t+MM (p, i) ) ; 

z=z+l ; 

end 

end 

return 
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o\o  o\o  o\o  o\o  o\o  o\o  o\o  o\o 


function [fill] =filtre (xf,  siglan)  ; 


Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
At  Air  Force  Institute  of  Technology 
Dayton-OH  USA  2002 


filtered  received  signal  calculator 
load  parameter; 

fiil=[] ; 

for  i=0:c-l; 

I=space3 (xf (l  +  i*subi :  (l  +  i)*subi,  :)  '); 
f ilt=I*siglan ( l+i*Fs : (l+i)*Fs) 
fiil= [fill  filt ' ] ; 

end 

return 
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o\o  o\o  o\o  o\o  o\o  o\o  o\o  o\o 


function [I] =space3 (M) ; 


Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
At  Air  Force  Institute  of  Technology 
Dayton-OH  USA  2002 


filter  weights  calculator 


aa=size (M) ; 

if  aa (1, 1) <aa  (1, 2) 

fprintf('you  know  that  spaceS  projects  into  columnsspace  of  M, 
right  ? ' ) ; 
end 

P=M*inv (M' *M) *M'  ; 

I=eye (length (M) ) -P; 
save  spaceSde; 
return 
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function [axx] =ambfun (x) ; 
load  sembold; 


g, 

0 

g, 

0 


%Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
%At  Air  Force  Institute  of  Technology 
%Dayton-OH  USA  2002 

g _ 

0 - 

^ ================================================= 

%detection  surface  calculator 


fdomain=fft ( [x  zeros (1, Fs*c) ]); %fdomain  contains  the  fft  of  the  signal 

vv=[0:20:400] ; 
dopp=length (vv) ; 


delta=totalt /Fs ; 

t= [ 0 : delta : totalt-delta]  ; 

axx=zeros (dopp,  2*Fs*c)  ; 
y=zeros (length (freq) , length (t) )  ; 

for  b=l:dopp; 

beta  =  1  +  2*vv (b) 73/10^8; 
factor  =  2*vv (b) *fc/3/10^8; 
dopp_signal= [ ] ; 
for  p=l:c; 
pt=t+ (p-1) *totalt; 

for  i=l : length ( freq) ; 

y(i,  : )=cos (2*pi* (freq(i) *beta+ factor) *pt+  (pi/2) *rn  (p, i) ) ; 
end 

dopp_sigp=sum (y,  1 )  ; 

dopp_sigp=dopp_sigp/max (dopp_sigp)  ; 
dopp_signal= [dopp_signal  dopp_sigp] ; 
end 

dopp_fdomain=f ft ( [zeros (1, Fs*c)  dopp_signal] ) ; 
axx(b, :)  =  abs (if ft (fdomain . *con j (dopp_fdomain) ) ) ; 

end 

save  ambfunde;% 
figure; 

me sh (lin space (-totalt*c, tot alt *c,2*c*Fs) *1000, vv, axx) ; 

xlabel ( ' Time  [ms] ' ) ; 

ylabel ( 'Velocity  [m/s]  '  )  ; 

zlabel (' Normalised  Ambiguity  Function'); 

return 


figure; 
clf ; 

orta=floor (dopp/2  +  1)  ; 
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if  vv(l)==0 
orta=l ; 

end 

newt=linspace (-totalt*c, totalt*c,  length (fdomain) )  ; 
range=150 ; 

sampdur=f loor (range*2/300000/ (newt (2) -newt (1) ) +1) ; 
if  sampdur>c*Fs 
sampdur=c*Fs ; 

end 

tstart=length (newt ) /2 ; 
tstop=  length (newt ) /2+sampdur; 

plot (newt (tstart:tstop) *1000, 10*logl0 (axx (orta,tstart:tstop) ) ) ; 
xlabel ( ' time  in  ms'); 

ylabel('the  mag  of  ambiguity  funcion  for  fd=0 ' )  ; 

title (' signal  auto  correlation  for  a  given  range  (here  150km)  for 
zerodoppler ' ) ; 


indis=f ind (xx (orta,  : ) ==max (axx (orta,  : ) ) )  ; 

figure; 

clf ; 

plot (vv, 10*logl0 (axx ( :  ,  indis )  )  )  ; 
xlabel  (' relative  velocity  in  m/s'); 

ylabel('the  value  of  ambiguity  function  for  zerotime  delay  along  the 
doppler  frequencies'); 

title ('this  plot  shows  basically  how  much  only-f req-delayed  versions  of 
the  signal  look  like  each  other'); 
save  cicimde 
return 
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function [ ] =gurultu (x)  ; 


g _ 

0 - 

g _ 

0 - 

%Written  by  Abdulkadir  Guner,  1st  Lt . ,  Turkish  Air  Force 
%At  Air  Force  Institute  of  Technology 
%Dayton-OH  USA  2002 

g _ 

0 - 

g _ 

0 - 

%gaussian  noise  calculator 

randn (' state ' , sum (100*clock) ) ; 
nn=randn (size (x) ) ; 

fprintf('SNR  for  normalizedcos/nn  will  be  -3  \n'); 
snr=input ( ' how  much  db  SNR  increase  do  you  want  in  SNR  ?'); 
if  isempty (snr) ==1 
ratio=l ; 

else 

ratio=10^ (snr/10) ; 

end 

nn=nn/sqrt (ratio) ; 
save  gurultude; 
return 
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